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The Influence of the Condenser on Microscopic Resolution 


By H. H. HOPKINS ann P. M. BARHAM 
Imperial College, London 


MS. received 13th April 1950 


ABSTRACT. An expression is found for the distribution of light in the image of two 
small apertures which are illuminated by a condenser of numerical aperture equal to s times 
that of the imaging objective. The resolving power of the system is written p=KA/(N.A.), 
and K is found as a function of s._ It is shown that precisely the same result obtains for 
both critical and Kohler illumination. 


oi INTRODUCTION 

HE theory of the resolving power of a telescope is generally based on a 
consideration of the image of a double star. Optically this is represented 
as the image of two closely spaced incoherent point sources of the same 
intensity. Each source produces a diffraction pattern in the focal plane of the 
telescope objective. ‘The image of the double star is then obtained as the sum 
of the separate intensities in the two diffraction patterns. The Rayleigh criterion 
of resolution requires that the minimum of intensity midway between the 
geometrical images is of the order of 20° smaller than the intensity at the maxima 
of the summed intensities. Conventionally the criterion is regarded as being 
satisfied when the geometrical image of one source falls on the first dark ring of 

the diffraction image of the other. 

In the case of a slit shaped aperture the intensity distribution in the image 
of a point source has the form J={sinx/x}*. ‘This gives a decrease in intensity 
of 18-9°% at the point midway between the geometric images. If a circular 
aperture is employed, the intensity distribution has the form J={2J,(z)/z}*. 
The decrease in intensity is then found to be 26-4°%. If « is the angular semi- 
aperture, and WN a refractive index, the two sources will be ‘resolved’ when their 
geometrical images are separated by a distance p, where 


Slit shaped aperture: = OUAN SUING ele te 0 set errr (1) 
Circular aperture: PN OLAWN Sioeie ts 5 Foes (2) 


It is significant that the smaller resolvable separation (1) is associated with a 
smaller drop in intensity. 

Disturbances from two distinct point sources are incoherent. In the case 
of a telescope imaging double stars, it therefore suffices to add the intensities 
of the separate diffraction patterns. This will not generally be true for the images 
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of, say, two illuminated pinholes in an opaque screen, for there some degree of — 
phase correlation will obtain. In general, therefore, it will be necessary to take 
phase relationships into account when adding the disturbances in the image — 
plane of a microscope. 

Martin (1931) considered the problem of the image of two small illuminated © 
apertures. His treatment was only approximate, but showed very clearly the | 
physical processes involved. ‘The case of critical illumination was considered. 
Firstly assuming a point source for the illuminant, and then extending the 
treatment qualitatively to consider an extended source imaged on the object | 
by means of a perfectly corrected condenser of aperture equal to that of the | 
objective, it was concluded that “the ultimate effect is likely to be similar to that | 
characteristic of two independent elementary sources”, the two apertures being © 
at a separation given by (2). In the present paper it is shown that this is, in fact, 
exactly true for both critical and Kohler illumination. ‘This accords with a_ 
result predicted by Zernike (1938). 

In our case an expression is obtained for the disturbance in the final 1 image | 
plane associated with the image of two pinholes illuminated by a condenser of © 
any given aperture when used in conjunction with a large source of uniform 
intensity. Circular apertures are assumed for both condenser and objective. 
The resolvable separation p is then written . 


p=KA/Nsin« 


instead of (2). K is found as a function of s, the ratio (N.A. condenser)/(N.A. 
objective). N.A. denotes, as usual, the numerical aperture equal to Nsin«. 

Recently a criterion of resolution due to Sparrow (1916) has been employed 
(Ramsay, Cleveland and Koppius 1941, Frangon 1950). If J is the intensity | 
at the point midway between the two geometric images, and w is a coordinate — 
measured along the line joining them, the images are said to be resolved when | 
o*I/dx?=0. It has been claimed that this gives a true ‘limit’ of resolution, 3 
since the intensity curve is then ‘flat’. However, the criterion is not unique, — 
since the equation 07//dx?=0 will frequently have an infinite number of roots. | 
Moreover the ‘ resolution’ indicated by this criterion is meaningless in some cases, | 
notably in finding the limit of detection of the microscope. The diameter of ° 
that smallest particle that can be ‘seen’ is found to be zero. With increasing » 
contrast sensitivity of the screen on which the image is received, the size of the : 
smallest detectable particle does decrease indefinitely. In a sense, therefore, the + 
result obtained using the criterion 0?//dx? =0 is valid. Nevertheless it is of little 
practical value. 

Because of this we have adopted a Rayleigh criterion of resolution which for - 
circular apertures requires the decrease in intensity to be 26-4% of the maximum | 
intensity. ‘The analysis is based on a consideration of critical illumination. It: 
is then shown that the same result obtains for Kohler illumination. 


> 


§2. CRITICAL ILLUMINATION WITH AN EXTENDED SOURCE 


We have first to define an optical unit of length. Suppose a length p in the} 
object plane has an image of length p’. Then the magnification is given by| 
m=p'|p. Ifa, «’ are the angles which a ray makes with the axis in the object 
and image spaces respectively, and these latter have refractive indices equal to 4 


i 
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N and N’, the magnification is also given by m=Nsina/N’sina’. Equation 
of these two expressions shows that the one coordinate length 

o= 7 pNsina == p'N’ sina’ bist (4) 
can be used to define equally p orp’. This is convenient since it is the variable (4) 
which appears in the functions expressing the amplitude in the diffraction 
distributions. 

Consider now an opaque plane containing two small apertures, P,, P,. The 
sizes of P, and P, will be assumed to be equal and both small compared with the 
Rayleigh resolution of the objective O, by means of which the two points are 
imaged at P,’, P,’. Let a, «’ in (4) refer to the objective O, and let p= P,P», 
=P, 'P,’.- Then 

2a 
Pa 7 Po(N-A.)o he ee (5) 


is the optical separation of Py and P,. (N.A.)y=Nsin« is the numerical aperture 
of the objective. 


Figure 1. 


If an element do of the source & is imaged at P(z,¢é) by the condenser C, 
the amplitude associated with the image of do will be that in an Airy disc with 


! centre at P. Let (1 =P Py, pp =P. ~Lhen the amplitude at P; due to do may 
_ be written 


Zi 
J, {= px(N.A.Jo} J,(s7;) 


2a STy 


_ where (N.A.)cis the numerical aperture of the condenser, and s =(N.A.)o/(N.A.)o- 
_ The variable 7, is defined by 7, =(27/A)p,\(N.A.)o. It 1s the optical measure of p, 
_ for the objective. The amplitude produced by do at P, is given by 


fi, SN at OTST foo a Sea ae (7) 
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where 7, is similarly defined for py. ‘The condenser is assumed to be free from 


aberration and perfectly focused. 


The diffraction images of the points P,, P, will be Airy discs centred on Py’, | 
P,’ respectively. ‘The amplitude produced at P’ by disturbances originating at_ 


P,, P, will therefore be given by 
J (a1) + Us J1(22) 


gy Las 


where 2), 2 are the optical distances P,’P’, P,’P’. The disturbances w,, uy are | 
coherent, since they derive from the same element do of the source. ‘To find the. 


intensity at P’ due to the whole source requires the integration of u? over the 
domain &. ‘That is 


Nees = | : {4 Oa + ty ae) do, i re (8) 


Expanding the expression in brackets in (8), 


rensi= {f urdo} XE) 4 Lf gtdo] EE 421 f ayuda} Hed 


B1 2 


Using the expressions (6) and (7) for uw, and u, the integrals in (9) may be evaluated. . 
To do this the domain & is assumed to be of infinite extent. Physically this }/ 
requires that the geometrical image of the source shall be large compared with : 


the distance P,P. 


Let P,, P, be at the points (—y,0), (+y,0) where 2y=z,. (2,¢) are the} 


coordinates of P. Then 
7,7 =y" + 2%+ 2yz cos ¢, 7.2 =y" + 2% —2yzcosd 
and the integral in the first term of (9) is 


20 2 1/ 
| atde= | ade ca ie re bag 


s(y? + 27 + 2yz cos )"/? 


The domain of integration has been formally extended to the infinite plane. . 


4 
} 


The integral remains unchanged if the origin is displaced to the point (—y, 0) 


It is then given by 
fe ee (sz ey ae aS 
ib ae ar { {> an 0 od 


eS) 


since | 
0 t 


21 3 
~ 
— 
EZ 
— 


bole 
° 


By a similar procedure the integral in the second term of (9) is found to be 


ie tido=nj[s. =) | ae (11) 


‘There remains the integral in the third term. 
This third integral is 


* susdo= ("eda (” Z1t0? +2" + 2y2 cos $)'%} Safs(y2-+.2?— 2ya cos $)1} 
Oo ae ha | ~ g(42 a o2'. Dye cowed onal -ale( nhac Seen 
Js 0 0 (+224 2yz cos $) s(y? +22 — 2yz cos $)! 
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which may be written 


; 7 22 J (sz) J,{s(3,2 +22 — 22,2 cos "1 
Gre . Ae ents ee eee ae ae 
| D3 pate | 0 ae eo 9 “Se $(89" + 3" — 2292 cos )!” “~ 


if the origin is displaced to the point (—y, 0), and 2y is replaced by 2». 

The integrand in (12) may be expanded with the aid of Neumann’s addition 
theorem for the Bessel functions (Watson 1942). Thus, consider the 
differentiation 


go o(R)} = — WY) absing, 
where R={a?+b?—2abcos¢}\?._ Neumann’s theorem has the form 


Jo(R) =Jo(a)F(b) +2. ¥ cos pd J,(a)J,(b). 


Differentiating this last with respect to ¢, and substituting for (0/0¢){J9(R)}, there 
is found the result 

aby die psin pd 

eo a ang FOI) 


If a=sz,, b=sz, substitution of this expression in (12) gives 


aia are = S psin pp 
| tata do = zs dz| ame I olsta)t ole) db. 
The integral in ¢ has the value 27 when p is odd, and is zero when p is even. 
re ice} J 68 Yf 
‘That is | UUs do = ue =p p50) | 1(sz)J (sz) oe 
= S” p=1,3... SZ 0 2 


The remaining integral in z is one of Lommel’s integrals (Watson 1942). It 
has the value zero for odd integral values of p, except p=1. In this case it is 
equal to 4. Hence 


: _ 1 2J,(s3) 
| tamedo= 5, atieltenearete (13) 
the remaining terms in the p summation being zero. 
Using (10), (11) and (13) in oe the intensity at the point P’ is found to be 


J Jy eu) 1°(2) | 7 24 also) J4(21) Ji(2) (14) 


an at S& By Be 


[(%, %y)= 


from which the common factor 7z/s? has been omitted. 

Some interesting conclusions follow at once from (14). If sq is a root of 
J,(s%) =0, other than sz) =0, the product term is absent. ‘I'he intensity I(3,, 2.) 
is then expressed by the sum of the separate intensities associated with P, and Py. 
P,, P, are therefore incoherently illuminated. In particular, this will be so, 
if s=1, when 2, is a non-zero root of J,(z%))=0: that is, when the points P,, P, 
are separated by a distance equal to the radius of any dark ring of the Airy pattern 
associated with the objective aperture. When P,, P, are on the limit of 
resolution (2), these points will therefore be incoherently illuminated if the 
numerical aperture of the condenser is equal to that of the objective. Under 
these conditions the resolution of two illuminated small apertures is given by 
the expression (2). 
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If the aperture of the condenser is made limitingly small, (14) becomes 


I(21) 22) = {ae 4 ae ee (15) 


1 oy 
since Lt {2J,(sz9)}/(s%o)=1. The points P,, P, are then coherently illuminated } 
s—>0 


no matter what their separation, 2). The significance of this in practice is best ¢ 
seen by noting that 2./,(s%9)/sz%)=0:880 when szg9=1. If s=4, points in a i 
circle of radius %)=10 will be very nearly coherently illuminated. ‘The diameter i 
of this coherently illuminated area (2%) =20) is approximately five times the limit 
resolution (%) =3-83) of the objective for two incoherent point sources. 

To investigate the intensity along the line P,’P,’, a coordinate x measured 


from the mid-point of P,'P,’ may be used. ‘The intensity at P’(x) is then given by y 


_ I(«+y) oe ee) 4.2 27 alszo) Ji (x+y) J,(x—y) 
(xy)? (a=) 5% Bay ey 


I(x) 


where y = 42, as before. 

It is a simple matter, using (16), to find the separations of P,P, which give: 
a decrease in intensity of 26-4% at O, when s has different values. In this manner | 
the value of K has been found for the range s=0-0—2-0. The curve of Figure 2} 


S = RATIO OF NAS 


Figure 2. Figure 4. 


shows the variation of K withs. If the influence of the aperture of the condenser + 
on resolution is ignored, K =0-61 for all values of s. This is the broken straight | 
line in Figure 2. The broken curved line shows K as a function of s on the basis : 
of the rule attributed to Abbé. According to this the effective aperture is the: 
mean of the apertures of the condenser and the objective. It leads to an absurd | 
result if the numerical aperture of the objective tends to zero. On the other hand | 
it gives a rough approximation to the present result if 0-5<s<1-0. When one 
remembers the other factors (such as scattered light, contrast in the object) | 
which influence the resolution of the microscope, it is not surprising that the Abbé / 
rule has been acceptable in practice. 
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§3. KOHLER ILLLUMINATION 


In critical illumination coherent light from an element of the illuminant is 
found in the object plane as an Airy distribution of amplitude. Thus a single 
element of the illuminant illuminates any two points in the object co-phasally, 
but with different amplitudes. In Kéhler illumination coherent light from an 
element of the illuminant illuminates two points in the object with light of the 
same amplitude but having a phase difference. The arrangement is shown 
diagrammatically in Figure 3. 


Figure 3. 


The source & is imaged by a condenser C, in the focal plane of a second 
condenser C,, which is further arranged to image a diaphragm D in the plane 
of the object S. Light from an element P of = therefore results in a plane wave 
falling obliquely on the axial region of the object plane S. ‘The object is 
illuminated by a set of plane waves whose normals lie within the cone subtended 
by C, at S. Let this cone be of semi-angle 8, and write NsinB=s(N.A.)o, 
(N.A.), being the numerical aperture of the objective, and N the refractive index 
of the space between C, and 5S. 

Consider two points in the object 5, such as P,, P, in Figure 1. ‘The plane 
wave falling on S whose normal makes an angle @ with the axis in the azimuth ¢ 
will have a phase retardation kNsin@.pcos¢ at P, relative to the phase at P,. 
‘The disturbances at P, and P, are thus given by 


Br =A, 
U, =A exp(—ikNsin@.pcos¢), | 


where A is the real amplitude at these two points. 

The condensers are assumed to be free from aberration. ‘They therefore 
satisfy the sine condition, and in consequence the rear unit surface of C, will 
be a sphere with its centre at the axial point of S._ If Ay is the amplitude of the 
incident wave at the unit surface H (Figure 4) of C,, and A is the amplitude of 
the emergent wave at the other unit surface H’, then the conservation of energy 
requires that A,*dhy=A*dh where dh, dh are elements of area. Now 


dh =dh,/cos 0, hence 
A tA COSN UF My we 8 ise ea (18) 


and the expressions (17) become 
ig, cos Ye, | 
Uy = A, cos!?Gexp(—ikNsin@.pcos¢). | 
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Let P,’, P,’ be the images of P,, P, (as in Figure 1) and let P’ bea point (2, 35). . 
Then the intensity at P’ is given by 


1%, 2) = [| u 


as in (8), but here w, vw, have the values given in equation (19). d@ is an element - 
of solid angle equal to sin@d@d¢, and the limits of integration are 0—2m for ¢ | 
and 0—f for 6. Further, since u, is complex, it is the squared modulus which 
appears in the integrand of (20). Hence corresponding to (9), (20) gives 


ee , Af 2 af 2 
I(%;, 23) = {f ip sind cos do dp {ja es ae 


B p2n : 
 Jal@r) Sil22) | | cos (kN sin 0. p cos g) . sin 6'cos 6 dO dd 
Oo’ 0 


By Bo 


J (21) a ee ie Soom (20) | 
2, 5 P| 


a constant factor A,? being omitted. 


The integral in the first term of (21) has the value wsin?8. The integral in | 


the second term is equal to 
B 
Os | sin 0 cos 0J9(kNp sin @) dO; 
0 


writing «x =kNpsin#@, and sz,=kNpsinf, this integral becomes 


27 sin? B /% ogee. 2I;(8%o)_ 
es {= J (x) de =n sin? B— 
Omitting the factor 7 sin? f, (21) is now 
I(%1) 2) = ats al eulee 7 24i(S%0) Fi(%i) Sul@) (22) 
Ber $25 = ce 


which is identical with ae Hence integrating the phases in Kohler illumination 


leads to precisely the same result as integrating the amplitudes in critical © 


illumination. 
This result is of considerable interest in practical microscopy. It shows 


that, given a broad source of uniform intensity, the appearance of the image will © 


be independent of whether Kohler or critical illumination is employed. The 


Kohler system has, however, two distinct advantages in other respects. There } 
is no difficulty in obtaining uniform illumination of the object, and the area of | 
the object in which light falls may be controlled by the diaphragm D. By this © 


means scattered light can be reduced with considerable improvement in the 
contrast of the image. 

It is unfortunate that critical and Kohler illumination are designated in French 
by the terms ‘éclairage incohérente’ and ‘éclairage cohérente’ respectively. 
The identity of the expressions (14) and (22) shows that the coherence between 


points in the object plane in the two cases must be exactly the same when the | 


condensers have the same numerical aperture. 

| 
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ABSTRACT. A direct photoelectric method allowed the reflectivity R and the trans- 
missivity T of partially reflecting films to be measured to an accuracy of 0:3°% of the incident 
light intensity. For a large number of silver films, deposited on glass by evaporation in 
vacuo, the values of R and T were measured for visible light of five different wavelengths, in the 
range of high reflectivities (R>0-7). The change of R as a function of time was also 
studied. With only moderate control of the conditions of surface cleaning and evaporation, 
the values of R were nearly always found to lie close to a ‘ standard’ curve giving R as a 
function of 7, the wavelength, and the age of the film. Possible causes of scatter from this 
standard curve were investigated. T'was found not to obey an exponential law as a function 
of the thickness. 

On the basis of these results, resolving power and brightness of etalon fringes were 
plotted as functions of JT. The conclusion is reached that, in most practical cases, the 
resolving power of etalons for visible light is limited by imperfections of the surfaces rather 
than by the reflectivity of the silver films, if their thickness is properly chosen. 


s1= INTRODUCTION 


N instruments based on the interference of multiple beams, transparent silver 
films of high reflectivity are widely used. The Fabry—Perot etalon and interfer- 
ence filters are among the most important examples. ‘The resolving power of 

these instruments is essentially proportional to 1/(1 — R), where R is defined as the 
ratio of the intensities of reflected and incident light. In increasing the value of R 
by increasing the thickness of the film, a practical limit is set by intensity require- 
ments, since the transmission factor T decreases rapidly. Beyond a certain limit, 
the resolving power can only be increased at considerable sacrifice of intensity. 
The accurate knowledge of the value of R which can be obtained for a given value of 
T at a given wavelength is therefore of practical value in interferometry. 

Though a number of observers have carried out direct measurements of 
reflectivities of silver films, our knowledge of the values of R and T in the range 
of high reflectivity (R>0-75) is far from adequate. 

The discrepancies between results of different observers and the scatter of the 
values obtained by any one observer are very large, and it is unknown how far these 
differences are caused by errors of the measurements, which are generally as great 
as 1%, and often more. 

Measurements of widths of interference fringes have given lower limits of 
reflectivities but, especially for high values of R, the results are considerably 
affected by imperfections of the plate surfaces. 

In order to study not only the values of the reflectivity as a function of trans- 
mission and wavelength, but also the scatter of the values of R for films deposited 
under similar conditions and the change of R with ageing, it was necessary to 
make the errors of measurement considerably smaller than these variations. ‘This 
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aim was achieved in the experiments described below. They were mainly 
intended to provide an answer to the following questions : (i) What value of R can 
be obtained for any given value of T between 0 and 0-2 and any given wavelength of 
visible light, if a good standard technique of surface cleaning and evaporation is_ 
used, and what is the amount of scatter about this ‘standard’ value of R ? (11) How, 
does the value of R change in the course of time, if the filmis kept under the same 
conditions as in ordinary spectroscopic use ? 

Regarding the influence of the numerous factors of surface conditions, vacuum 
and speed of evaporation, some useful information was obtained, but a compre- 
hensive study of this question was not attempted. 


Sal HE OP rr eC Aree Ae OD 


The principle of the method employed is simple and well known, but the: 
accuracy of the results depends on some details which are not quite obvious and 
will therefore be described (Figure 1). 

The plate P, carrying the film whose reflectivity is to be measured, is mounted l 
on a spectrometer table, with the axis of rotation passing through the reflecting 5 
surface. ‘The metal plate holding P has a circular aperture Q of diameter 1-1 cm. . 
The light source S is in a fixed position, 200 cm. from P. A photo-multiplier M is 5 


20 cm. ao 200 cm. > 


Figure 1. ‘The optical system. 


mounted on the movable arm of the spectrometer. In position M,, the photo- « 
multiplier receives light passing from S through Q directly when P is removed, or ° 
after passage through the plate P when this is inserted. In position My, it: 


receives light reflected from P at a small angle of incidence (about 15°). This 


light can be regarded as coming from the virtual image S’ of the source. The : 


distance of S’ from M,, is the same as that of S from M,. The intensities observed 
in the three experiments are therefore in the ratio 1: 7: R, provided the following 
conditions are fulfilled: (a) The light source and the aperture of the photocell 
must be so small that any ray passing from one to the other intersects the plane of P 


well inside the stop Q. (6) The light source must emit uniformly in all directions | 
within the solid angle subtended by the aperture of the receiving system, and the } 
latter must have an isotropic response. These conditions are especially important | 
since S’ is the mirror image of S, so that lack of uniformity in the angular / 
distribution of the emitted light will have different effects in the two positions of M. | 


(c) Light scattered by the mount of the reflecting plate, by the reflecting surface 
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itself, and by any other parts of the apparatus must be negligible. (d) The 
reflecting plate P must be sufficiently plane. (e) All geometrical adjustments have 
to be sufficiently accurate: the plane of P must pass through the axis of rotation, 
and the plane defined by S, M, and M,; must be normal to the axis of rotation, 
and the lines SM, and SM,,; must pass approximately through the centre of Q. 

In order to fulfil conditions (a) and (6), light from a 24-watt bulb B, connected 
to a large battery, was focused on a ground glass G which was fitted with a circular 
stop of diameter5mm. A similar ground glass G’ and stop were fixed in front of 
the photo-multiplier. 

‘The surface of the brass plate holding P was covered with black velvet paper. 
When M was set only slightly off the angle of reflection, the intensity fell to zero, 
showing practical absence of scattered light including light scattered from P itself. 
The intensity measured as a function of small horizontal or vertical displacements of 
M from the reflection position showed a symmetrical curve with a flat top whenever 
all adjustments had been carefully made. 

The photo-multiplier, type R.C.A. I.P.21, was connected directly to a galvano- 
meter. By varying the distance of the light source from the photocell and applying 
the inverse square law, it was found that the galvanometer deflection was not strictly 
proportional to the intensity, and all readings had to be corrected accordingly. 

The effective wavelength range in each experiment was a function of the 
colour filter used, of the spectral intensity distribution of the source and of the 
response curve of the photosensitive layer. The latter has a maximum at about 
4100 a. and the sensitivity decreases rapidly in the near infra-red. From spectro- 
scopic examination of the filters, the effective wavelength was estimated. 

The Table gives a list of the filters used, the estimated limits of the 
effective wavelength bands, and the approximate centre of gravity of each band. 


Ilford red 6400-7200 a. 6800 a. 
Ilford green 5100-5400 a. 5200 a. 
Ilford blue-green 4700-5200 A. 4900 a. 
Wratten No. 48 4400-5000 a. 4700 a. 
Wratten No. 35 3900-4500 a. 4200 a. 


The violet, blue and blue-green filters were used in conjunction with a cell 
containing copper sulphate solution. 

The colour filters F (Figure 1) were placed very close to the small aperture in 
front of the photo-multiplier. When they were placed further away, the small 
local variations of density of the filters were found to cause errors up to 0:5 %,. 

When a measurement was repeated, after removal and re-setting of P and 
readjustment of the instrument, the value of R was generally reproduced to within 
0-3°%, and disagreements greater than 0-5 °/, were practically never found. ‘The 
estimated limits of error of a single measurement of R was thus about + 0-003. 
With the precautions described, all possible systematic errors should have been 
practically negligible, and the stated limits of error are believed to represent the 
absolute accuracy of the values of R. ‘The accuracy of measurement of T is a 
little greater, so that the value of R+T for each film should rarely be in error by 
more than +0-005. 

In the method described, the value of R was measured at an angle of incidence of 
15°, whereas in the Fabry—Perot interferometer the angle of incidence 1s practically 
zero. For opaque metal films, the change of R with the angle of incidence can be 
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calculated from the values of the optical constants and is found to be of opposite _ 


sign for light polarized in the plane of incidence and at right angles to it. For the 
average of the two planes of polarization, the difference between the values of R for _ 
0° and 15° is quite negligible (of the order of 0-001). For transparent films, the 
theoretical formulae are complicated, but for the comparatively low transmission 
used in our experiments, the dependence of R on the angle cannot differ much from _ 
that for thick films. This was indirectly confirmed by some measurements with 
plane-polarized light; the differences of the values of R for the different planes of 
polarization were small and close to the values predicted by the theory for thick 
films. 

The error made in regarding the values of R measured at an angle of 15° as the 
reflectivity for normal incidence is below 0-002 for dense films, but it may arise to 
about 0-01 for the lowest reflectivities measured. 


§3. THE PREPARATION OF THE FLIEMS | 

The evaporation chamber consisted of a large brass tank evacuated by means of a | 
mercury diffusion pump of a rated speed of 25 litre/sec., but the liquid air trap and 
the duct to the chamber may have reduced the speed to about half this value. The 
pressure in the tank was measured by means of a Penning gauge. It was between 
3 and 5 x 10->mm. Hg during the evaporation. ‘These values refer to calibration 
with air. 

Two tungsten filaments of thickness 0:7 mm., each carrying two beads of silver, 
could be heated independently. A few turns of thin platinum wire were wound 
round the V-shaped parts of the filaments, according to common practice, in order 
to cause the silver to adhere to the filaments. When a fairly good vacuum had been 
reached, air was admitted to a pressure of a few millimetres of mercury, and a 
discharge from a spark coil was passed through the tank for 5 to 10 minutes before 
the pumping was resumed. ‘This treatment was repeated once or twice. ‘The 
actual evaporation lasted about 5 minutes. 

In each evaporation, three sets of two glass plates were used, at distances of 18, 
20,.and 22cm. respectively from the filaments. This method allowed the 
comparison of films deposited under exactly the same conditions of evaporation. 

The glass plates were of good optical quality, accurately plane and parallel, 
and of size 4-1cm. x 3:2cm. Some of the plates were polished lightly with rouge 
before they were taken into use, so that the surfaces were similar to those of new 
plates, but this process was not repeated after the plates had once been covered with 
silver. ‘The plates were washed in nitric acid, then rinsed with tap water and 
dried with cotton wool of good commercial quality. ‘The surface was considered 
satisfactory when breathing produced a uniform film. 

For the ageing tests, the plates were kept in a glass container, protected from. 
dust by a loosely fitting cover. ‘The container was left in a Laboratory room which 
was free from tobacco smoke, flames and chemical fumes. 


§4. THE MEASURED VALUES OF R AND 7 


In the final set of experiments, the results of which are reported here, 18 


glass plates were coated with silver in the way described and each plate was. 


measured four times, first immediately after removal from the evaporation tank, 
then after intervals of 3, 8, and 22 days from the time of deposition. Each time 
the values of T and R were measured for five different wavelengths. Six plates. 
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were left to be re-measured after 68 days, while the other 12 plates were cleaned, 
re-silvered and measured again when new and after 22 days. 

For the measurements of these 30 films immediately after evaporation and 
after 22 days, the values of R+ T are plotted in Figure 2. Within the range of T 
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Figure 2. Measured values of R+ 7 as function of 7. 


plotted, these graphs contain all the measurements, without any selection. Only 
one film gave values so much below the curves (by about 5%) that they were not 
plotted for technical reasons, while for another film, for the same reasons, only the 
measurement in the red for the new film was included. 
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Since 7 changes with the wavelength, the same films do not appear at the same | 


! 


groups of films can be seen to form similar ‘constellations’ of points in the graphs | 


values of the abscissae for different curves. If this fact is taken into account, some 


of different wavelengths. The group of 6 points which appears at 6800 a. between 


T =0:06 and T=0-09 and at 5200 a. between T=0:12 and T=0-16 forms an — 
example of this fact which shows that the scatter of the points in the graphs is | 


mainly due to genuine differences of the properties of the films, and only to a small 
extent to errors of measurement. 

Most points lie quite close to a curve showing a decrease of R+T with 
increasing transmission, i.e. decreasing thickness of the film. ‘The absorption 
A =1-R-T is thus found to increase with decreasing film thickness. Except for 
the violet where the scatter is somewhat larger, most points lie less than }°% off 
the curve. 


The curves of Figure 2 were used in plotting R as a function of T in Figures | 


3 (a) and 3 (bd). 
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Figure 3. R as function of T. 
(a) for new films, (6) for aged films (22 days). 


The loss of reflectivity of the films as a result of ageing depends on the wave- 


length of the light and on the thickness of the film. Examination of Figure 2 | 


and of other data shows that for dense films R decreases less rapidly than for thin 
films. In order to obtain a simplified picture of the process of ageing, on the 
basis of statistical evidence, the relative values of R+ T were averaged for the 


six films which were re-measured after 68 days, i.e. the average of (R+T)/(R+T) 9 


for these six plates was plotted as a function of time in Figure 4, where (R + T)p 
is the value of R+ T for the new film. Since T' changes little and is much smaller 


than R, this plot is almost identical with a plot of R/Ry. The six films were of | 
different thicknesses, but all of them, expressed in terms of interferometry, were | 


fairly dense films for red, and fairly light films for violet light. 
The values for red and green show that the rate of decrease of R becomes small 
after three weeks. ‘The points for 4200 a. cannot be fitted well by a simple curve; 
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examination of the individual results for the six plates and for others which were 
measured after 50 days makes it appear likely that the ageing curve for violet light 
is, in fact, not simple. In any case, it is certain that the reflectivity decreases with 
age more rapidly for shorter wavelengths. 


1-00 | | 
68004. 
099 i | | 
& S 
: 0-98 x 
S 
N 
Hon 
ona 4200A. 
0-96 
0-95 | | | “ 
0 10 20 30 40 50 60 70 
Days 


Figure 4. Relative change of R and R+T with age. 


A number of measurements made outside this systematic research confirms the 
conclusion that, in general, the rate of loss of reflectivity. decreases with time ; one 
pair of etalon plates, whose properties had not been measured accurately 
immediately after deposition of the films, were found to have values of R+ T of 


4000 5000 6000 7000 
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Figure 5. Optical density D=log1/T as function of wavelength, 


‘ 0-97 and 0-98 for red light after having been in use for over 6 months. In some 
» films, the reflectivity was actually found to have increased after a few weeks, but in 
». these cases, the initial values of R + T had been abnormally low. 

| Figure 5 shows the optical density D =log 1/T plotted as function of the wave- 
‘length. The points forming each curve are averages of several films of very 
i similar density. 
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§5. DISCUSSION OF THE RESULTS 


Before this systematic research was undertaken, the instrument described in 
§2 had been used* for measuring the values of R and T for a considerable number | 
of silver films deposited on etalon plates, most of which were made of fused silica. | 
Most of these values agree well with the curves in Figure 2. 

Since the scatter of the points in Figure 2 is much greater than the errors of 
the optical measurements and thus represents genuine differences 1n the properties 
of the films, it appeared desirable to establish the causes of these differences. 
But the large number of variables which might possibly influence the properties of 
the films allowed only tentative conclusions to be reached. 

(i) Influence of conditions of evaporation. If vacuum conditions and speed 
of evaporation, within the range of variation in these experiments, had an important 
influence on the properties of the films, any two films deposited simultaneously at 
the same distance would tend to be more similar to one another than to other films | 
of similar density. No effect of this kind was apparent. However, in comparing | 
films of approximately equal density, but deposited at different distances from the | 
target, some evidence was found for the smaller distance to produce slightly higher - 
reflectivity. This effect, if real, could be interpreted as due to different speeds of | 
evaporation.t / 

(ii) Influence of the type of glass. Of the 18 glass plates used, some consisted 
of glass appearing distinctly green when viewed end-on, while others appeared 
quite colourless. ‘The results obtained with these two kinds of glass did not show 
any systematic differences. Silica plates, as mentioned before, also gave similar 
results. 

(iii) Influence of surface conditions. Of six pairs of glass plates, each pair | 
was coated together at the same distance, and the properties of the films were : 
measured several times. The films were then removed with acid and the plates : 
re-coated in the same way ; in one case this procedure was repeated a third time. . 
In six of these seven repeats, the plate which gave the higher value of R, compared | 
with the ‘standard’ curve of Figure 2, in the first coating, also gave the higher ° 
value in the subsequent coating. ‘The differences in question were of the order of [ 
0-01. It thus appears that some surface conditions which are not effected by the : 
cleaning process, are often responsible for producing a ‘ better’ or ‘ worse’ ’ 
coating. In a few cases, imperfect cleaning was, no doubt, the cause of too low ; 
values of R, and on rare occasions even of much too low values. 

Previously published direct measurements of reflectivities of transparent films : 
(Romanowa, Robzow and Pokrowsky 1934, Goos 1936, Krautkramer 1938, , 
Strong and Dibble 1940) show so much more scatter than the values reported here : 
that only rough comparison is possible. 

The work of Strong and Dibble (1940) covers mainly the range of low and | 
medium reflectivity and their measurements of R are only relative and hardly | 
accurate enough for application to interferometry. The existence of two different | 
types of silver films below a certain thickness, as reported by these authors, may be : 
connected with our observation that thinner films are less stable than thicker | 


* These measurements were carried out in conjunction with Mr. L. C. Bradley and Dr. G. W. | 
Series, and the results were reported at the International Conference on Optical Properties of | 
Thin Films in Marseilles in April 1949. 

} Note added in proof. R.S. Sennett and G. D. Scott (J. Opt. Soc. Amer., 1950, 40, 203) report | 


a marked increase of reflectivity with increase of speed of evaporation. 
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ones. Our measurements do not extend to sufficiently thin films to allow direct 
comparison. 

The measurements by Krautkramer (1938) fall into the range of higher reflecti- 
vities, but were made with a small number of films only. Within their limits of 
error, stated as + 0-01, they are compatible with our results. 

Among others, Tolansky (1946) and Dufour (1948) have given values of 
reflectivities, but lack of details does not allow comparison. 

The reflectivity of opaque silver films, made by evaporation, was measured by 
Edwards and Petersen (1936). ‘Their values for yellow and green light (R =0-985) 
are in very good agreement with our results in the limit of T=0 (see Figure 2). 

Values of the reflectivity of silver films were derived from measurements of 
the width of interference fringes produced with etalons of very small spacings 
(Bright, Jackson and Kuhn 1949). ‘This indirect method gave lower limits of 
R+T which are perfectly consistent with our direct measurements. The increase 
of R+ T with increasing thickness of the films was, in these indirect measurements, 
disguised by the imperfections of the etalon plates, an effect which becomes more 
important as R increases. 

For opaque films, the reflectivity can be calculated from the optical constants, 
but the results of such calculations have generally been found to be too low. If 
the values from the latest measurements of the optical constants of freshly deposited 
silver (Hass 1946) are compared with our limits of R for T=0, this discrepancy 
still exists. ‘The measured values of R exceed the derived values by amounts 
increasing from 0-02 in the red to 0-06 in the violet. 

For transparent films, the formulae of the theory are rather complicated. 
Calculations by Barns and Czerny (1931) show a rate of decrease of R+ T with 
increasing T which is similar to that found in our experiments, but their absolute 
values are consistently lower. 

Though the thickness of the films was not measured, the results plotted in 
Figure 5 show that the intensity of the transmitted light does not vary exponentially 
with the thickness of the film ; if the exponential law of absorption were valid, 
the ratio of the optical densities for any two films would be the same for different 
wavelengths. Figure 5 shows that these ratios are quite different for red and 
violet. Some indication of this fact was found by Bright, Jackson and Kuhn (1949), 

One can demonstrate this effect, which is, of course, theoretically not unex- 
pected, by a simple qualitative experiment; a pile of glass plates, each carrying a 
thin silver film, shows a deep violet colour in transmission, very different from the 
blue of a single, dense film. 


§6, APPLICATION TO INTERFEROMETRY 


The resolving power of the Fabry—Perot etalon can be defined 1n various ways, 
and a number of authors have tabulated it as a function of R.A very convenient 
definition (Bright, Jackson and Kuhn 1949) considers two lines as just resolved 
if their distance is equal to the instrumental half-value width; the resulting 
intensity curve shows a depression of about 20° for lines of equal intensity, as in 
the Rayleigh condition for gratings and prisms. A simple expression for the 
resolving power can then be derived if in Airy’s formula the sine of the phase 
angle is replaced by the angle itself, near the maxima of the intensity. ‘This approxi- 
mation is accurate enough, as long as the reflectivity is not much below 0:8. If 
the resolving power is expressed, in the usual way, as the product of the order of 
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interference n and a factor N which can be considered as the effective number of 
inter fering beams, the latter is found to be 


N=nvV/Ri(1—R). 


With the use of the measurements described above, the value of NV is plotted 
as a function of 7, for new and for 22 days old silver films, in Figures 6 (a) and 6 (0). 
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Figure 6. Effective number of beams N and intensity J of etalon fringes. 
(a) for new films, (6) for aged films (22 days). 


In the same Figures, the intensity of the maxima, J ={T/(1—R)}’, is plotted. These | 
graphs make it possible to find the value of 7, for a given colour, which gives the | 
best compromise of the requirements of resolving power and intensity. If, for 
example, an intensity of the fringes of 0-3 of that without any etalon is required, fora 
several weeks old film, Figure 6 (6) shows that for red light 1/110 of one order can be 
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resolved, and for green light 1/40 of one order, provided the etalon plates are 
perfect. ‘The required transmission factors are 0-015 and 0-04 respectively, for 
these two colours. If the value of T is to be measured with light of any other 
colour, the curves in Figure 5 allow the necessary conversion to be made. 

It must be remembered that the formula J={7/(1—R)}* gives the correct 
intensity only for monochromatic light. If the spectral width of the line under 
investigation is comparable with the instrumental width, the intensity is lower. 

The surfaces of commercially obtainable etalon plates are generally not good 
enough to allow a smaller fraction than 1/50 of one order to be resolved, and less 
than this for plates of large size. Our measurements show that, certainly for 
wavelengths of over 5000 a., a considerable increase in resolving power could be 
achieved with ordinary silver films, if etalon plates with more perfect surfaces were 
available. 
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ABSTRACT. Recent interest in lead sulphide, selenide and telluride as photoconductors 
makes information on their absorption spectra valuable. It is found that the absorption 
spectrum of each material is characterized by strong absorption in the visible and ultra- 
violet region, with a further band in the infra-red. The effects of temperature and oxidation 
are also studied. 

In PbS the forbidden energy region between the filled and empty bands is found to be 
about 1:3 ev. in width, not 0:4 ev. as previously suggested. Hence photoconductivity in 
PbS cannot be associated with the main lattice. The corresponding energy gaps in PbSe 
and PbTe aré found to be about 1:05 ev. and 0-9 ev. respectively. 


Si INDGRODUCTION 
HE absorption spectra of insulators and semiconductors are usually 
characterized by very strong absorption at short wavelengths (usually in 
the ultra-violet region) with a sharp edge on the long wavelength side. 
This absorption is generally assumed to be due to electronic transitions between 
the highest filled energy band and the lowest empty energy band of the solid. 
The absorption coefficient is usually of the order of 10° or 10® cm™}. 

The lattice absorption edge is often followed on the long wavelength side by 
an absorption band or tail of much lower absorption coefficient, of the order of 
10 cm-!. This absorption is generally dependent on the conditions of sample 
preparation and is thought to be due to impurities, cracks and other defects in 
the crystal lattice. 

The absorption spectra of PbS, PbSe and PbTe are not exceptions to these 
general conditions, though the long wavelength tail absorption bands have 
considerably larger absorption coefficients and bandwidths than those usually 
encountered in sulphides. 

The photoconductive properties of PbS and its associated compounds are 
well known (Sosnowski et al. 1947, Moss 1949, Chasmar 1948). The spectral 
response of all three materials is characterized by a sharp cut-off at long wave- 
lengths occurring at about 3 for PbS and about 5y for PbSe and PbTe. In 
all cases the long wavelength limit moves to longer wavelengths on cooling. 
Early theoretical work (Sosnowski et al. 1947) suggested that the photo- 
conductivity of PbS was due to absorption in the bulk lattice. This implies 
that the minimum energy difference between the highest filled band and the 
conduction band of PbS is about 0-4 ev., and less for the selenide and telluride. 

The photoconductivity of PbS and its associated compounds is largely 
determined by suitable treatment with oxygen (Sosnowski et al. 1947, Schwarz 
1948). Recent work in this laboratory (Gibson and Moss 1950) together with 
other unpublished work by Mr. A. S. Young, has shown that oxygen treatment 
has very marked effects on the photosensitivity spectrum. Thus in the case of 
Pb'Te photoconductivity is limited to wavelengths less than 34 without oxygen 
treatment, moving to 5p on oxidation. This and similar work suggests that 
photosensitivity in these materials is not primarily a property of the bulk lattice. 


| 
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§2. EXPERIMENTAL METHODS 

Lead sulphide and lead selenide layers may be prepared by two methods: 
evaporation 7m vacuo or chemical deposition from colloidal solution (Kicinski 
1948). Layers of PbTe have not, as yet, been successfully deposited chemically 
and hence only evaporated PbTe layers have been examined. Layers were 
mounted on discs of polished quartz, artificial sapphire or KRS5 (a mixed crystal 
of thallium bromo-iodide) depending on the wavelength range under examination. 
Chemical methods of deposition gave layers of the greater reproducibility and 
uniformity of thickness. Within the limits of chemical analysis (1°%), PbS 
layers deposited chemically are pure PbS. 

The chief errors in absorption measurements on thin films are as follows: 
Firstly the layers consist of small microcrystals (of the order of 1000. in size) 
and transmission between grains may be more important than transmission 
through grains. ‘This effect is easily detected as an apparently constant absorption 
with wavelength and may be eliminated by careful layer preparation. In 
particular, chemical layers are superior in this respect to evaporated layers. 
A second important source of error is due to scattering of light by individual 
grains of the layer. ‘This effect may be eliminated by collecting the transmitted 
radiation in a spherical mirror of high aperture. 

The third remaining error is due to reflection. As the refractive index of 
PbS, PbSe and PbTe is very high (about 4:0 for PbS and over 5-0 for PbTe) 
the reflection coefficient is high and reflection represents a most serious source 
of error. When the absorption within the specimen is large (say 99% or greater) 
the error due to reflection can be ignored. When the percentage absorption 
falls to low values, however, on the long wavelength side of the lattice absorption 
edge, reflection introduces an apparent absorption which may be greater than the 
true absorption. 

To eliminate reflection errors two specimen layers were prepared 
simultaneously under exactly similar conditions, except that one layer was 
arranged to be thicker than the other. The transmission through each layer 
was measured concurrently at each wavelength. The chief reflection, which 
occurs at the air—specimen interface, is the same for both layers. Errors due 
to reflections at the specimen—mount interface and mount-air interface are only 
equal in each layer if the absorption in the specimen is zero, but this 1s where 
reflection errors would otherwise be most serious. When the layers differ in 
transmission by a factor of five, the remaining reflection error is about 1%. 

If the two layers have thicknesses d, and d, respectively, the absorption 
coefficient being ~, then 


J,=J,(1—R) exp (—«ad,) and J,=J,(1—R) exp (— ad) 
where J, and //, are the intensities of the transmitted beams, J, the intensity of 
the incident beam and R the reflection coefficient of the air—-specimen interface. 


In J,/J5 
Hence — pS 


In the presentation of the results log/J,/.J, has been plotted. ‘To obtain the 
absolute value of « the difference in layer thickness must be known. ‘This is 
obtained by weighing, but high accuracy is not possible. Values of « are 
indicated alongside those of log./,/-/5. 
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The errors inherent in this method that are not accounted for are as follows: 

(i) Interference between reflected and incident beams. This effect would 
lead to an apparent absorption varying cyclically with wavelength. In fact no 
such effect was observed. This was thought to be due to the corrugated nature 
of the specimen surfaces. 

(ii) If both the mount and specimen absorb at any particular wavelength, _ 
then the amount absorbed by each mount is not equal. This error can be 
arranged to be negligible. 

(iii) In spite of identical preparation conditions a pair of layers may not, in 
fact, be identical. Tests showed that this source of error was unimportant in 
chemical layers, but occasionally evaporated specimens had to be abandoned 
because of this feature. 

Another method of eliminating reflection errors was also examined. ‘This 
consisted of measuring the transmission through a ‘black body’ type cone of 
Pyrex glass coated with PbS and collecting the transmitted radiation in a 
spherical mirror of large aperture. ‘The absorption curve obtained over the 
limited wavelength range available, (0-3 u to 3-41) was the same as that obtained 
by the difference method. ‘The cone method was not very suitable for evaporated 
layers and was later abandoned. 

The absorption measurements were made using a double prism monochromator 
with interchangeable quartz, LiF and NaCl prisms, using a thermocouple, a PbTe 
cell or a photomultiplier as detector. A wavelength resolution of 0-05 » or better 
was generally employed. 

To study the absorption spectra at low temperatures both layers of each pair 
were totally immersed in liquid nitrogen in an unsilvered quartz dewar vessel. 
For some measurements layers were deposited directly on to the inner surface 
of the dewar vessel. ‘The accuracy of these measurements at low temperatures 
was not high, say 20%, due to the constant boiling of the coolant and water 
condensation. 


EXPERIMENTAL RESULTS 


93) THE ABSORPTION SPECTR UM Oi LE AUS Ulla Eis iE, 


3 (1) General Results 


The absorption spectrum in the range 0-2 to 7-0 of a typical chemically 
deposited layer of PbS is given in Figure 1. If care is taken to use a standard 
deposition method the layers are quite reproducible. Deviations from the 
standard method cause small changes in the absorption spectrum, particularly 
around 0-9. Layers may be prepared with a very marked dip at this wavelength 
which separates the lattice absorption and the tail into two distinct bands. 

Evaporated layers of PbS give very similar absorption spectra to those of 
chemical layers. ‘The chief differences occur (a) in the long wavelength tail 
band which may differ in magnitude by as much as a factor of 3 and (b) in the 
wavelength region around 0-6. The original PbS available for the preparation 
of evaporated layers is, however, of doubtful purity. Thus two samples of | 
‘pure PbS’ available commercially have been found to contain 45° and 55% 
of PbSO, respectively. The chemically deposited layers probably represent 
the purest specimens available. 
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The chief characteristics of the absorption spectrum of PbS are easily 
distinguished, the lattice edge being very marked. The absorption coefficient 
in this region is about 210° cm, which is a reasonable value for any 
semiconductor (Mott and Gurney 1940). The most surprising feature of the 
absorption spectrum is the marked tail band which extends beyond 6u. The 
absorption coefficient in this band is about 10* cm~! and is therefore much larger 
than that found in semiconductors. The significance of this result will be 
discussed later. 


log (A/a) 


x 
o 
wa 


Absorption Coefficient (cm=! 


02 03 05 0:7 10 i) 20 3 5 7 10 
Wavelength (microns) 


Figure 1. The absorption spectrum of chemically deposited lead sulphide. 


3(u) The Effect of Oxygen Impurity 

The importance of treatment with oxygen to obtain photosensitive layers of 
PbS is well known. 

Oxygen may be introduced either (a) during the preparation of the original 
material, or (b) during evaporation, or (c) after the layer has been deposited. 
Generally speaking only processes (b) and (c) have marked effects on the 
photosensitivity. 

Chemical layers of PbS are normally baked in air for about 30 minutes at 
about 150°c. in order to produce photosensitivity. ‘This corresponds to 
oxidation process (c) above. Initial measurements indicated that such treatment 
made no significant change in the absorption spectrum. Bakes up to 16 hours 
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duration were tried. This result was not altogether unexpected as von Hippel 
et al. (1946) have reported the same result for thallium sulphide, another material _ 
normally sensitized by baking in air or oxygen. | 

A significant effect due to oxygen was observed when care was taken to reduce — 
all extraneous effects to a minimum. In particular the layers were not removed — 
from the monochromator mounting for baking. Evaporated layers were prepared _ 
in vacuo and their transmission measured at each wavelength before air at 
atmospheric pressure was admitted by breaking a small glass seal. 

The change in absorption observed was very small (about 3°% or less) but — 
almost doubled on baking in air. The magnitude of the change was wavelength 
dependent, being most marked at the wavelengths 1:0 and 2-3, in the tail 
band. In the tail absorption band the absorption decreased on admission of 
oxygen. | 

In view of the small change in absorption to be observed the errors in this — 
measurement are large, but the reality of the effect is beyond doubt. The most _ 
surprising feature is that the absorption coefficient in the tail absorption band © 
decreases on oxidation, when an increase would be expected if new impurity bands _ 
are being introduced. It may be significant that previous workers (Sosnowski 
et al. 1947) have noticed peaks in the photoconductivity spectrum of PbS at 
about 1 and 24 which they ascribe to oxygen treatment. ‘These wavelengths 
agree fairly well with the figures given above. 

It has been stated by Wilman (1948) and others that as much as 10% of 
PbO.PbSO, may be formed on the oxidation of PbS films. It is almost certain 
that the main absorption band of this material lies in the far ultra-violet outside 
the range of measurement. 


3 (ii) The Effect of Temperature 

The absorption spectrum of PbS at room temperature and at 77°. is shown 
in Figure 2, the abscissa being the reciprocal of the wavelength. The most 
obvious effect of cooling is to shift the lattice absorption edge to longer wavelengths. 
A number of determinations of the magnitude of the shift have been made, the mean 
value indicating a shift of about 6 x 10-4ev/°kK. This compares with the measured _ 
shift in the spectral limit for photoconductivity of 4-8 x 10-*ev/°K., and is in the 
same direction (Moss 1949). 

A lattice edge shift to longer wavelengths on cooling is not usual in semi- 
conductors. ‘There are, in general, two factors which determine the shift of the 
lattice edge with temperature. As the lattice contracts on cooling ions move 
closer together, interaction increases and allowed energy bands broaden. In 
addition, however, electron interaction with the thermal vibrations of the lattice 
decreases on cooling and allowed energy bands narrow (Radkowsky 1948). It 
is not possible without detailed calculation to determine which effect will be the 
more important in PbS. ‘The energy change due to contraction is very approxi- 
mately proportional to energy bandwidth, whereas the shift due to lattice vibration 
effects is independent of bandwidth. As the former effect appears to be the more 
important in PbS it is concluded that the allowed energy bands in PbS are very 
wide. In general this would be expected im a sulphide of high dielectric constant. 

The long wavelength tail absorption band of PbS also undergoes a change 
on cooling. The band appears to narrow appreciably, at least on the long 
wavelength side, and the peak of the band moves to longer wavelengths. The 
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total area under the tail band is not appreciably altered. It should be pointed 
out that measurements at 77°K. could not be made at wavelengths greater than 
4-2 « because of absorption in the quartz dewar vessel. In addition measurements 
around 2-7 were impossible due to the condensation of water on the cooled 
surfaces of the dewar vessel. ‘Thus the accuracy in the 3 « and 4p region is not 
very high. 
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Figure 2. The absorption spectrum of lead sulphide. 
A at 77°x., B at 290° x. 


3 (iv) Other Measurements on PbS 


(a) In view of the lack of marked effects due to the exposure of evaporated 


~ PbS films to the air, layers were prepared in vacuo on KRS5 discs and subse- 
+ quently opened to the air. As KRS5 is transparent to very long wavelengths 
~ (to about 40) it was possible to examine the absorption of PbS out to 16 (the 

' limit set by the monochromator and detector). No major absorption bands in 


the region 7 to 16m were observed. The measurements set an upper limit 


to the ‘ Reststrahlfrequenz’ in PbS. 


(b) Attempts have been made to measure the absorption spectra of single 
crystals of natural galena. ‘These experiments were not very successful as it 


was impossible to obtain crystal specimens sufficiently thin and yet free from 


holes. [The measurements obtained indicated a similar structure to that shown 
in Figure 1, the lattice edge being easily observable. 
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Due to the above difficulties it was not possible to observe changes in 
absorption with crystal direction by transmission measurements. As an 
alternative method crystals were split and polished and mounted on a table 
rotatable through 270°. The reflectivity of the specimens was then measured 
with plane-polarized infra-red radiation at various crystal angles. No significant 
changes in reflectivity with crystal direction were observed. 


$4. THE ABSORPTION SPECTRUM OF LEAD SEVENTDEE 


The absorption spectrum of a chemically deposited layer of this material 
is given in Figure 3 where the absorption coefficient is plotted against the 
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Figure 3. The absorption spectrum of chemically deposited lead selenide. 


logarithm of the wavelength. ‘The spectrum is similar to that of PbS except | 
that it is shifted bodily to longer wavelengths. The shift is about 0-4 ev. in the: 
region of lattice absorption. 

The lattice edge of PbSe, like that of PbS, moves to longer wavelengths | 
on cooling. ‘he shift is about the same as, or rather less than, that of PbS. | 
The expansion coefficients of PbS and PbSe are almost identical but the allowed | 
energy bands are likely to be somewhat narrower in PbSe. According to Moss | 


(1949) the shift in spectral limit of photoconductivity of PbSe is exactly the same 
as that of PbS. 
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yo Loa eABSORP LION SPECTRUM OF LEAD TELLURIDE 

This material cannot be obtained by chemical deposition and hence only 
“fevaporated films may be studied. PbTe can, however, be prepared by heating 
the elements in sealed and evacuated quartz tubes and hence a product of high 
purity (particularly as regards oxygen) obtained. ‘This represents an important 
“jadvantage over PbS and PbSe. 

The absorption spectrum of a layer of PbTe prepared by the evaporation 
in vacuo of a powder sample manufactured in vacuo is given in Figure 4. The 
‘general features are similar to those of the sulphide and selenide. 
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Figure 4. The absorption spectrum of lead telluride (material prepared in high vacuum). 


The effect of oxygen on the absorption spectrum of this material can be studied 
both in the bulk material and in films. As for PbS, the absorption in the long 
i | wavelength tail band decreases slightly when air is admitted to vacuum-prepared 
‘films but the effect is smaller than that in the sulphide and too small for adequate 
measurement. 

To study the effect of oxygen on the bulk material, specimens of PbTe have 
_ been prepared by heating the elements in low pressures of air and in open crucibles. 
The effect of such treatment is illustrated in Figures 5 and 6. ‘The absorption 
‘at short wavelengths is not appreciably affected (the main peak at 0-55 1 remains 
‘entirely unchanged in position or magnitude) but marked changes occur at 
| longer wavelengths. A new absorption peak appears at about 0-9, and the 
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tail absorption decreases in magnitude on its short wavelength side. The latte 
effect results in an apparent shift to longer wavelengths of the tail band maximu 
and the minimum between the main bands. 

This comparatively violent oxidation (the material examined in Figure 6 
was heated to over 1,000°c. in air) is clearly similar, except in magnitude, 4 
oxidation of sulphide and telluride films by exposure to the air. Similar thoug 
less marked effects can be obtained by evaporation of unoxidized PbTe sample 
in low pressures (of the order of 10-2 mm.Hg) of air. It seems likely that al 
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Figure 5. The absorption spectrum of lead telluride (material prepared in low pressure of air). 


the more marked differences in absorption spectra between PbS, PbSe or PbTe 
samples of different origin are due to oxidation effects. 

Very intense oxidation of PbTe results finally in complete decomposition 
and the formation of a grey-white product. By analogy with PbS, the material 
is thought to be PbO.PbTeO,. The only marked absorption of this material 
commences about 28004. and extends below the limit of measurement. It 
seems likely that the absorption bands of PbO.PbSO, would be at even shorter 
wavelengths. 

The shift in the lattice absorption edge in PbTe on cooling is even less than 
that of the selenide and its magnitude proportionately uncertain. A reasonable 
value would appear to be a 1:5 x10-*ev/°k. According to Moss the shift in 
photoconductive limit is the same as for the sulphide and selenide but recent work 
in this laboratory by Mr. A. S. Young has shown that any shift between 2 and 
7 x 10-*ev/°K. may be obtained in PbTe, depending on sample preparation. 


The Absorption Spectra of PbS, PbSe and PbTe 765 


yo" THE ABSORPTION SPECTRA OF PbS, PbSe AND PbTe: 
GENERAL DISCUSSION 


6(i) Lattice Absorption 


The lattice absorption bands of each material almost certainly correspond to 
lectronic transitions between the highest filled energy band and the lowest 
pty band, that is, the transfer of an electron from S, Se or Te ions to Pb ions. 
he long wavelength part of this band may correspond to exciton transitions 
Mott and Gurney 1940), but there is no evidence in favour of this view. 
Extrapolation indicates that the lattice absorption coefficient of PbS will 
each zero at about 0-94, corresponding to a minimum energy of about 1-3 or 
4+ev. ‘This value is considerably greater than the value of 0-4ev. suggested 
y Sosnowski et al. (1947) and suggests that these authors were mistaken in 
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i Figure 6. ‘The absorption spectrum of lead telluride (material prepared in air). 


associating photoconductivity with lattice absorption. It would appear that the 
photoconductive spectrum is in fact determined by impurities or other effects. 

It is of interest to note that the sensitization spectrum of optically sensitized 
PbS layers (Gibson 1949) is almost exactly coincident with the main lattice 
absorption. Previously sensitization was thought to take place at specks of free 
sulphur but the similarity between the sensitization spectrum and the photo- 
conductivity of free sulphur is not very close. It is now clear that a more likely 
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explanation of optical sensitization involves absorption in the main lattice and thi 
subsequent filling of electron traps lying just below the conduction band of thi 
solid. 

The threshold energy for lattice absorption falls progressively from : 


sulphide to the telluride, as might be expected. The magnitude of th 
absorption coefficient increases through the series. The latter effect mus 
be associated with decreasing energy bandwidths as the spectra move to longe} 
wavelengths. Due to this progression almost all the main absorption band oj 
PbTe lies within the limits of measurement and an estimate of the ee 
may be made. The half width of the absorption band is about 4 ev., which mus 
be approximately equal to the sum of the widths of the two allowed energy band 
involved. : 


\ 


6 (ii) The Long Wavelength Tail Absorption Band 


| 


By analogy with other photoconductors and semiconductors this band may 
be ascribed to impurities in interstitial or substitutional lattice positions. Theré 
are, however, a number of reasons for rejecting this explanation, as follows. 

(a) The absorption coefficient is too large, being about 10* cm™ in PbS ana 
largerin PbTe. Absorption bands of the magnitude and width observed would 
correspond to an impurity concentration of the order of 107° per cm*. PbS! 
PbSe and PbTe are highly polarizable materials and therefore it may be argue 
that they are capable of taking up very large amounts of impurity without seriout 
lattice distortion. Notwithstanding this, it seems unlikely that a chemically 
deposited film of PbS, which is never heated, should contain about 100 times a: 
many impurities as are usually found in other sulphides. 

(b) If photoconductivity is associated with impurities and the excitation oy 
electrons from impurity levels, the absorption bands observed cannot correspond 
to such transitions, as the absorption spectra and photoconductivity spectra ar 
not the same. At room temperature absorption in PbS extends to 6 or 7p 
photoconductivity to 3p. 

An alternative explanation is to ascribe the observed absorption band te 
electrons in the conduction band of the solid, the transition involved being t» 
the next highest empty band, which probably overlaps the conduction band 
‘The measurements of Eisenmann (1940) and others suggest that the concentrationi 
of electrons in the conduction band of PbS is of the order of 1019 or 10° at room 
temperature. ‘This number appears to be largely invariant with temperature 
(Dunaer 1947). Thus, PbS, PbSe and PbTe are semi-metallic (Chasmar 1948) 
These electrons would give rise to an absorption band at long wavelengths of abou: 
the magnitude observed. ‘The area under the band would be largcly unalteree 
on cooling, though the energy distribution of the electrons may alter. Ne 
photoconductivity would arise from this type of transition and hence the lack 0: 
photoconductivity at 6 is accounted for. ‘The effect of oxygen in diminishing 
the absorption band may be interpreted as the capture of electrons by oxygel 
atoms-on the surface of the crystallites. 

If this explanation is correct, photoconductivity cannot arise from an increast 
in the number of electrons in the conduction band as the efficiency would be to¢ 
low. Schwarz (1949) has suggested that photoconductivity is a surface phenomen: 
arising from adsorbed films of oxygen. A theory of this type has been considere¢ 
for some time in this laboratory. Each crystallite of PbS (n-type) is assumed t 
have oxygen atoms adsorbed on its surface which capture electrons from th 
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anterior of the crystallite and produce a negative surface charge. Illumination 
fmay eject electrons from the surface oxygen states, the surface charge would be 
reduced and conduction facilitated. ‘There is increasing experimental evidence 
in favour of this view but conclusive evidence is still lacking. 
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ABSTRACT. An analysis has been made of a large number of hot cathode low voltage 
arc-type discharges in mercury vapour and hydrogen at low pressure by the cold probe 
methods of Langmuir and Mott-Smith. It has been verified that the space potential in 
the plasma near the filament is positive relative to anode potential, and it has been shown 
that consistent data for the space potential and concentrations of the slow (ultimate) electrons 
are obtainable from analyses of probe characteristics in both retarding and accelerating | 
fields for electrons. The way in which the fast (primary) electrons develop a random) 
motion has been traced by progressively confining them within the main discharge, firstly 
by using a cylindrical gauze anode with open ends, secondly by using a solid anode with : 
open ends and lastly by using a solid anode with end caps. 


SiS SIN LRODUE TON 


accelerated through a double space-charge sheath and produce a highly | 

conducting plasma further out. The collision processes occurring in the é 
plasma have been investigated in detail by Langmuir and Jones (1928), who have 
studied how the primary electrons have their directions altered and give up energy | 
in various types of collisions with the gas molecules. They were able to distinguish 
three main groups of electrons in the plasma; the primary electrons from the : 
filament, the secondary electrons ejected from atoms and produced by loss of | 
energy of the primary electrons, and the slow ‘ ultimate’ electrons which are < 
mainly responsible for the high conductivity. The tube used had a straight ‘ 
filament cathode mounted axially in a closed cylindrical metal box, with the end : 
caps insulated from the curved wall. By varying the gas pressure, Langmuir and - 
Jones made an investigation* of the transition from the velocity distribution - 
principally studied, which can be approximately described by the three groups, to ° 
the velocity distribution which is encountered under many other conditions : 
(Langmuir and Mott-Smith 1924), in which a single group is present, often witha 
Maxwellian distribution of velocities. "he present author has found that some - 
simple modifications of the electrode system used by Langmuir and Jones permit | 
the demonstration of the transition in another way. ‘The modifications consist in A 
the use of different types of anode, which progressively confine the faster electrons } 


IE an hot cathode arc discharge at low pressure electrons from the cathode are 


more within the general body of the arc plasma. This shows rather directly how } 
the overall velocity distribution tends to become Maxwellian, although both the | 
atomic process involved and the current voltage characteristics of the collectors | 
are too complex to permit full theoretical analysis. 


* See particularly Figure 13 of their paper. 
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§2. APPARATUS AND EXPERIMENTAL METHODS 
Most of the experiments were performed in mercury vapour saturated at a 
/temperature of approximately 20° c. A few experiments performed with hydrogen 
1 at the same pressure gave similar results. No allowance has been made for’ 
thermal effusion. The tubes had an axial tungsten filament cathode 0-2 mm. in 
‘jdiameter, used as a bright emitter. The length of the filament between its 
supporting leads varied from 6-0 to 6-5 cm. as replacements of it were made. The 
‘fall of potential along the filament was approximately 8 volts. The anode of the 
‘first tube (Figure 1) was made of nickel gauze (with 25 wires/cm., each wire 0-15 mm. 
if diameter) of length 4+-2cm. and diameter 2-3cm. In the second tube the nickel 
/) gauze was replaced by nickel sheet of the same dimensions, and in the third tube 
jthe nickel sheet was furnished with two flat end caps through which the filament 
jand probe leads passed. Two cylindrical molybdenum probes were used. The 
first F, was situated at about the middle of the filament, and the second F, at one 
)end in such a position that its exposed length always lay within the perpendicular 
»from the corresponding end of the anode on the filament. Both probes were 


0-1mm. in diameter, F, being 0-8cm. long and F, 0:4cm. long. They were 
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Figure 1. Diagram of peered and leads for tube 1. Gauze anode, AA; filament, f, held taut 
by sprung metal strip, S, at left-hand end; probes, F, and F;, with glass shields omitted. 
The four leads on the right passed to a common pinch. The lead to F, passed through a 
separate seal at the left. 


+ shielded by fine glass tubes except at the exposed ends. ‘The probes were parallel 
7 to the filament, and about 3mm. from it. They were put close to the filament in 
/ order to minimize complications from plasma oscillations, which do not usually 
+ become pronounced until farther out (Neill 1949). By reversing the direction of 
the current through the filament, information could be obtained about the 
‘discharge at four different regions along its length from the negative end. The 
electrodes were mounted axially in a glass tube 4-5 cm. in diameter and 15 cm. long, 
; with an appendix containing a drop of pure mercury. ‘The tube was kept 
a continuously evacuated by a diffusion pump. 

Most measurements were made with the tubes passing temperature-limited 
¥ currents from the whole of the filament. A few were made with the current not 
i quite saturated. No significant differences were found between the two sets of 
“results. With tubes 1 and 2 the plasma could be seen to bulge beyond the open 
» ends of the anode, in some instances filling the whole of the tube. With tube 3 
+ there was occasionally a small penetration of plasma past a filament lead. ‘The 
characteristic curves were obtained using standard circuits as described, for 
+ example, by Greeves and Johnston (1936); any runs in which variation of probe 
/ potential caused a detectable change in anode current were rejected. ‘The general 
' shape of the characteristic curves was that described by Langmuir and Jones. 
. For sufficiently negative voltazes of the probes a small positive ion current was 
- received wnich did not vary much with voltage. When the probe voltage was 
| PROC, PHYS, SOC. LXIII, IO—B 3G 
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raised to that of the neighbouring part of the filament a rapid change in current 
occurred, due to the reception of primary electrons. With increasing positive 
voltages more primary electrons were received. On approaching the anode 
‘potential a further rapid increase of electron current to the probe occurred, 
due to the reception of the ultimate electrons. In no case was it possible to’ 
distinguish clearly between primary and secondary groups. 


§3. ANALYSIS OF RESULTS 

As the probes were situated close to the filament it was essential to know that 
they were actually in the plasma. This was considered likely, from the appear- 
ance of the discharge and from the fact that the cathode sheath was probably much 
nearer the cathode (Druyvesteyn and Penning 1940). It was confirmed by the | 
analysis of the electron currents received by the probes when their potential was : 
near that of the anode. Since the ultimate electrons were present in much greater 
numbers than the primary electrons, the latter could be allowed for on this part 
of the probe characteristic by approximate graphical extrapolation. ‘The residual 
currents of ultimate electrons were analysed in two ways, first, by plotting the | 
currents on a semi-logarithmic scale against the probe voltage (i.e. (log7z, V) 
graph), and second, by plotting their squares against the probe voltage, on a 
linear scale (i.e. (22, V) graph). Both methods permit the determination of the 
space potential and of the electron concentration (Langmuir and Mott-Smith 
1924). Inthe second method it is necessary to apply two corrections: first, for - 
the position of the space potential, using a value for the mean electron energy 
obtained by the first method, and second, for the distortion of the linear : 
[(current)’, potential] graph through incomplete orbital motion near the space : 
potential. Reasonably consistent results were obtained by the two methods, as is : 
shown by some typical data given in the Table. In all cases the space was # 
positive with respect to the anode. 


Plasma Data 


Ultimate electron 


Tube Pressure Arc current Space potential (volts) concentration (10°/cm?) 
and probe (10-% mm.) (ma.) from (log7,V) from (z7,V) from (logi,V) from (2?,V) 
fy 1235 160 TS 7-4 25-6 20-9 
il, 1% 1235 176 70 7A 0:3 0:2 
De le 1:24 220 5:0 5:0) PUD 18-8 
A), lls 1:24 200 2°2 2-0 io 0:8 
Shy 14 1-04 347 4:5 455 18-6 16:8 
Sah 1-04 340 1:7 1°7 0-4 0-4 


The arc voltage drop was 25, and the filament current 4-85 amp. in each case. 


The distribution analysis outlined above suggested that if the primary electrons ; 
were progressively confined within the inter-electrode space by partial reflection : 
from the surface of suitably designed anodes, the overall distribution would tend to } 
become more nearly Maxwellian. Tubes 1, 2 and 3 (see §2) were designed and | 
used one after the other with this end in view. In tube 1 a number of the faster | 
electrons were being lost to the main body of the plasma because they passed out / 
through the wires of the gauze. In tube 2 there was no penetration of the anode by | 
fast electrons, although not all were reflected from the metal (Langmuir and Jones | 
1928). In tube 3, loss of electrons through the ends of the anode was prevented in | 
addition. Results of the type anticipated were obtained. 
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Consideration was given to the possibility that the space potential was near the 
athode potential, corresponding to a point on the (log 7, V) graph at which the 
ate of increase of primary electron current with increase of probe voltage begins to 
iminish rapidly. ‘This alternative interpretation of the (log 7, V) graph has been 
ejected from the close analogy between the discharges studied in this investigation, 
nd those studied by Langmuir and Jones (1928) and Compton and Eckart (1925) 
n which the space potential was near that of the anode. 

An attempt has been made to determine the velocity distribution amongst the 
rimary electrons from analysis of the (log 7, V) curves with the probe potential 
ot far from that of the cathode. Qualitatively the distribution, particularly for 
ube 1, clearly consists of a large component of drift motion with a superimposed 


190 
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}Figure 2. Semi-logarithmic plots of electron current 7 on arbitrary scale against probe voltage V 
; relative to anode. Curve (b) is displaced 12-5 volts to the right and curve (c) 22:5 volts 
to the right. The space potentials are at S. 

| 

|smaller random motion. Methods for analysing such curves have been developed 
by Langmuir and Mott-Smith (1926). ‘Three complete semi-logarithmic plots of 
‘electron current against probe voltage are shown in Figure 2, for tube 1 (curve a), 
‘tube 2 (curve b) and tube 3 (curvec). ‘The discharge conditions were: for curve a, 
pressure 1-35 x10-*mm., are current 160_ma. ; for curve 4, pressure 1-24 10-3 mm., 
are current 220 ma.; for curve c, pressure 1:04 10-3 mm., arc current 347 ma. In 
‘each case the arc voltage drop was 25-0, and the filament current 4-85 amp. 
‘Further data are included in the Table. ‘The curves show a transition from the 
form in which the primary electrons are well separated (a) towards a form in which 
the semi-logarithmic plot is more nearly straight over the greater part of its length 
negative to the space potential. ‘This type of transition has been found consis- 
tently in some 100 characteristic curves for F,; it occurs inall cases by simultaneous 
flattening of the left-hand hump H, and extension towards more negative voltages 
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of the approximately straight part corresponding to the ultimate electrons near the i 
space potential. Since the existence of a straight-line curve for (log 7, V) shows 


to have here an example of the degeneration of the separate groups towards a) 
single group having an overall Maxwellian distribution and equivalent temperature. : 
This is what would be expected through the progressive enclosure of the discharge, 
with diminishing leakage from the cylinder of primary electrons, and, possibly to a| 
less extent, of ultimate electrons. The characteristic curves have proved to be too\ 
complex to permit quantitative analysis, but it does appear that not only is there 
a modification of the primary group, but in addition, an extension towards higher} 
energies of the ultimate group, and that, could the discharge conditions be varied 
widely enough to follow these processes still further, a single group with a definite 
temperature (much above the gas temperature because of the presence of electric) 
fields) would be obtained. 
Confirmation of the essential correctness of this interpretation of the data fon 
F, is afforded by the results obtained with F,. When this was at the negative end) 
of the filament, the characteristic curves for tubes 1 and 2 were similar to one 
another, whilst with tube 3 only a slight straightening occurred. ‘The similarity o 
the characteristic curves for tubes | and 2 is due to the magnetic field of the filament 
current. ‘This bends the electron beam emerging from the filament sheath intos 
the cylinder, and therefore robs the region surrounding this end of a considerable! 
number of primary electrons. ‘This bending occurs with both tubes, and although 
the laws of electron reflection are unknown in detail and probably vary consider-4 
ably, the increased reflection from the solid anode in tube 2 can scarcely result iny 
the redirection of many electrons towards the negative end. Hence the overall! 
velocity distribution there remains practically unchanged. In tube 3, again, 
little alteration of the characteristic curves would be expected because the axial 
component of the velocity of the primary electrons is small. 
When F, was at the positive end of the filament the probe characteristics were 
found to change in type much as they did with F,, for all three tubes. Thisé 
again is to be expected, because at this end the magnetic field tends to make thew 
primary electrons pass out of the cylinder. : 


ae 
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“ABSTRACT. A method is described for the measurement of the viscosity of molten 


© metals at temperatures up to about 900° c. ‘The metal is contained in a cylindrical crucible 


? which is allowed to oscillate about its axis under the control of the torsion in a suspension 
“wire. The theory of the method is considered. 

The viscosity of zinc has been determined in the temperature range 430° c. to 480°c. 
)and the results are given. 


SV ieuN DRO DUC BTN 


HERE is not a great deal of information available on the coefficients of 
| | viscosity of molten metals. Still less is known of the viscosities of alloys 

and the effect, upon the viscosity of one metal, of the addition of small 
quantities of another. 

Most of the well-known methods for the determination of the viscosity of 
substances which are liquid at room temperature are unsuitable for measurements 
fon liquid metals. The use of the oscillating disc method, in which a solid of 
» revolution such as a disc, suspended by a wire, is allowed to oscillate about its 
"axis in the liquid, against the restoring torque provided by the torsion in the 
wire, has been described in some detail by Stott (1933) who applied it to the 
investigation of the viscosity of molten tin. ‘The chief disadvantage of this 
+ method is that contamination of the surface of the liquid, by oxide or otherwise, 
leads to a serious dragging effect on the rod connecting the disc to the suspension 
wire. It is difficult to eliminate this contamination in experiments on substances 
such as molten zinc. The method of enclosing the liquid inside a hollow sphere 
} attached to a bifilar suspension and allowed to execute torsional oscillations has 
| been developed and used with success by Andrade and Chiong (1936). 

They discuss the calculation of the coefficient of viscosity from the observed 
| damping of the oscillations and show that the method can be made to yield 
| accurate results. Andrade and Rotherham (1936) give an alternative method 
of determining the viscosity by cbserving the forces necessary to maintain the 
oscillations at a constant amplitude. 
. The choice of material for the construction of the hollow sphere presents 
problems when it has to contain molten metal at high temperatures. ‘The 
material must be immune to attack by the molten metal, and must be capable 
of being formed to the required shape without too much difficulty. Sintered 
alumina is a substance which withstands high temperatures and resists chemical 
attack, but hollow spheres which can be sealed with the metal inside are not 
easily made from this material. Accurately made cylindrical crucibles of sintered 
alumina are, however, readily available, and it seemed desirable that an 
investigation should be made of the possibility of determining the viscosity of 
molten metals by observing the oscillations of a cylindrical containing vessel 
Icft open at the top end. ‘This method has been suggested by Knappwost (1948). 
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For a cylinder it is more difficult than in the case of a sphere to obtain 
manageable expression from which the viscosity can be calculated in terms o 
the dimensions of the apparatus. It is possible, nevertheless, to use the method 
for absolute determinations. The theoretical considerations outlined belo 
indicate how this may be done, and also provide the basis for a calibration formula, 
which can be determined by using liquids of known viscosity. 


| 
| 
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The oscillation of a closed cylinder of length short in comparison with its: 
diameter has been considered by Meyer (1891), and viscosity determination 
made on the basis of his analysis by Mutzel (1891).* | We shall consider here th 
case of a long cylinder open at one end. Let the cylinder, containing liquid, 
execute rotational oscillations about its axis (which we take to be the z-axis, in th 
vertical direction). Let (wu, v, w) be the velocities of the liquid at the poin 
(x, y, 2) in the (x, y, 2) directions, and let p (x, y, z) be the pressure. Then, 1 
(X, Y, Z) are the components of the body forces and v the kinematic viscosity 
n/p where 7 is the coefficient of viscosity and p the density, the general equations; 
of viscous flow are 


Gus Ou en Cues 1 op ee 

ae Secs Kat win, 

dv dv dv OU lop 3 

GPa, ue ae a ee a aivethowatcs (1) 
Oe ee ge ee 
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to which must be added the continuity equation 
du dv dw 
a ale ay ai oe =. oe ee (2) 


Assuming that there is no motion in the z direction, that there are no body‘ 
forces, and that the velocities are small enough to enable us to neglect theé 
non-linear terms, the equations reduce to 


ou _—s— lop 


ae ee ee 2 
ay Ber ei at F 
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Ou dv 

stat | (4) 


Let us assume now that each particle of fluid moves in a circle with a velocity} 
depending only on the distance r from the axis, on g and on ¢, and put 


u=—py, V= x, | 

where ¢ is a function of 7, x and t. The tangential velocity at a distance r is then mp / 
and the radial velocity zero; % is in fact the angular velocity of the liquid. | 
* We are indebted to the referee for pointing out that work on the determination of viscosity } 


by observing the oscillations of a closed cylinder containing liquid has been done by Okaya and |) 
Hasegawa (1933) and Okaya (1936). | 
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+) Multiplying the first of equations (3) by y, the second by x and subtracting, we 


have, in cylindrical polar coordinates, when there is no azimuthal variation of yb 
ay 30h OM _ Lay 
Ge oes ay fie °) 


We suppose that the system has settled down so that the liquid has the same 
period and decay as the containing vessel. 


Then pb =dew 
and if a is the radius of the cylinder, c the height of the liquid in the cylinder, 


+ and the origin is taken at the centre of the base of the cylinder, 


w=De™ at 2=(0 and “at? =a, 
0 


ae =(0 at z=c. 


Equation (5) then becomes 


Oe AO On 30D. O°@ 
wee pies oe Byer (6) 
with the boundary conditions 
d= at z=0 and at Cle 
real fi 
be =( at z=c. | ”) 
dz 


a is complex; in fact «=—fB+v2y, where B=28/T and y=27/T, 5 being the 
logarithmic decrement and T the period of oscillation. The solution of (6) 
satisfying the boundary conditions (7) ts 
- a® J (mr) 2 20 y J ,(R,r)m? cosh L,(2 —c) p 

r Jima) r  Jo(kya)R,(R,? —m") cosh l,c’ 
where —m*=«/v and J, and J, are Bessel functions of order 0 and 1 respectively. 
The summation is taken over the positive roots k,, of J,(k,,a) =0, 1, being given 
by /,2=m?—k,?. The damping torque on the system due to the viscous drag 
of the liquid on the walls of the cylinder is 


kes ol ye Oe, | 
G=2n\a iE a), 2 I. =F) cada reece (9) 
Of _ nt Of _ tO 
where Erle as nd ae ae ae 
which, on evaluating with reference to (8), becomes 
‘- mca*J o(ma) Ae = 
G=2nn \- ~ J,(ma) + 2m QA £4 Lk, 2(k,2— me) i poogon (10) 


where =e is the angular velocity of the cylinder containing the liquid. 
We can write ‘VY =d6/dt and 6=e™/« where @ is the angle through which the 
cylinder has turned from some zero position. 

Equation (10) can be expressed as 


G=L— IIA 
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and the equation of motion of the oscillating system, assuming that the external 
damping due to causes other than the viscous drag is small, is 


a0 dé 
=e fe =H ee eee 1Z 
: dt? dt ees oe 


where J is the moment of inertia, and f the restoring torque per unit angle of 
rotation. 
As is shown by Andrade and Chiong (1936), (12) may be written as follows on 
putting L=M+iN and remembering that «= —f+1y. 
Ion? +(M+iN)a+f=0, 
I(p?—?)—-BM-yN+f=0, 
—2Ipy—BN+yM=0, 


i.e. m= 6 {pha tihy N=7{ et - ih. 


Since 8B =28/T, y=2n/T, f=47?I/T,2 where T and Ty, are the periods of the | 
full and empty cylinder respectively. 


Zid (alae 
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Pores (13) 
eee a 
TRIN ie 
where A=8/y =8/z7. 
Since A is small, the expressions for M and N become 
Zio 2al (-f 
M="r(75+1); N= "3 (a3-1). ee (14) 


These equations contain only measurable quantities, so that M and N are known. 
Simple considerations show that if the external damping is small but not 
negligible, it is approximately accounted for by subtracting the logarithmic 
decrement observed when the cylinder is empty from that observed when the 
cylinder contains liquid, to obtain the value 5 in the expression for M. The 
apparatus should, however, be designed to keep the external damping as low 
as possible. 

The value of the coefficient of viscosity 7 in terms of the logarithmic decrement 
6 can theoretically be obtained by combining equations (10), (11) and (14) to give 


Zod J (ma) 
ma +1) = Real t of eS ee 
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or 
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In order to use this, it is necessary to simplify the right-hand side. This we 
now proceed to do. 


Therefore m = (x) a —1) (1 a ey : 
ay 
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To a first approximation therefore 


a \ : : 
and ma=a( 7) (1—7)=p(1—2) say. 
Tv 
It is convenient to put ¢=(1—7) and observe that 


JS) 2  So(S) 
JQ) &  Sa(2) 


Since, in our experiments p is large, we can use the well-known asymptotic 
expansions for J,(f) and J,(¢), viz. 


110~(2)"{(-vbe) (6-8) d(-9) 
10~(2)"{(t+ ee) oo6-9) fee 9} 


Now cos (¢—7/4) =e“e“-7/9/2 and sin(¢—2/4) =e"e™"/21 on neglecting e™ 
compared with e”, so that 
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which reduces to 


gerAeiisre) again neglecting terms in ele compared with unity. 
8u—3 | 
AG (43) 
Hence Oo = +7 
: patel arena ea 
| J (5) J (5) 
and (a 
(5) J(5) 
8ut1 
Equation (15) may now be written 
PAR of aS Bi 8ut 1 } : | 
ai ine = 2 yy ee K ein ai 16 
TF (73 ae 1) 2arq| ca {2 (5) + K(p) (16) 
tanh Jc 


x as D 
where K(y)=Real part of 22 24 bA(ak,)?— ©] 


men 


If the height of the liquid in the cylinder is not small compared with the diameter 
of the cylinder, K() is much smaller than the other term on the right-hand side 
of (16). Physically, this means that most of the viscous drag takes place at the 


Tales F| 
Peet 
Cees { ? “ << = a ~ othe 
i Ny ree } ER! ea ae ee iad 


778 M. R. Hopkins and T. C. Toye 


curved surface of the cylinder as compared with the end. We can proceed now ~ 
by using (16) to find an approximate value of 7 and of course p, neglecting K(), 
and then using the approximate value of 4 to evaluate K(u) and so by a process — 
of successive approximation arrive at a value of 7 satisfying the equation. It is . 
found in practice that, after first obtaining an approximate value of 4, it does not 
involve much extra labour to compute J,(¢)/J:(¢) without neglecting terms — 
in 1/£2, and thus to obtain numerically a more accurate expression than that given 
by the term containing {2— (84+ 1)/(8—3)} in (16). 

Calculation of K(w) for various values of » shows that it is very nearly a 
linear function of ». It is worth tabulating this function if the apparatus is to | 
be used unchanged for a large number of determinations. 

An alternative procedure is to write (16) in the form 


ar? 9T 
8 (7s + 1) =AnT+ BunT+ re 


3771/2 
=AnT+ B(npT)?+C KE ] 


and to determine the constants A, B and C by observations on liquids of known 
viscosity. The cubic equation in 7 can be solved by the usual methods. 


$3. APPAR AWS 


The apparatus is described with reference to the sectional diagram of Figure 1. 

The silica furnace tube K, 3 ft. 6 in. long and 5 in. internal diameter, is wound 
at E. ‘The winding is arranged so as to ensure that a uniform temperature can 
be maintained over a length of the furnace safely in excess of that occupied by 
the sintered alumina crucible H which contains the molten metal under 
investigation. [he dimensions of the crucible are : height 10cm., internal 
diameter 3:9cm., wall thickness 1mm. Below the crucible are heat insulating 
bricks F which are supported by means of a long steel tube which rests on the 
lower steel endplate T of the furnace. The temperature of the furnace is 
controlled by a Cambridge regulator, and measured by means of a carefully 
calibrated platinum—platinum-rhodium thermocouple, used in conjunction with 
a potentiometer graduated in hundredths of a millivolt, but capable of being 
read by estimation to a tenth of this value. For the type of thermocouple used, 
one millivolt corresponds roughly to 100°c. The regulator is used to connect 
one or other of two alternative sources of current to the furnace winding. By 
careful adjustment of the relative magnitudes of the currents available from 
these two sources, and by maintaining the cold junction of the thermocouple 
at 0c, in a vacuum flask containing melting ice, it is found that the temperature 
can be controlled to +$°c. ‘This variation takes place over a period which is 
long compared with that required to make a determination of the viscosity. — 
The furnace tube is surrounded by the lagging J. 

The crucible H is supported in a light steel cradle attached to a steel rod | 
which passes through the steel tube supporting the bricks F in the upper half — 
of the furnace. At its upper end the steel rod carries an aluminium disc D, | 
having a graduated flange, which can be observed through a glass window in 
the brass cover G. The whole oscillating system consisting of the crucible, 
the connecting rod and aluminium disc is suspended by a phosphor bronze wire, 


| 
‘i 
i 
| 
4 
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ten inches long and one hundredth of an inch in diameter, fastened to a brass 
rod L which fits closely into a brass sleeve, the joint being lubricated with vacuum 
grease. ‘The suspended system is located vertically by the brass bush R clamped 
to the rod L. Slight rotation of R starts the oscillations. 

The upper part of the apparatus, consisting of the brass cover G, the upper 
endplate 'T of the furnace, the upper insulating bricks F and the whole of the 
suspension system, can be removed as one integral part, the only vacuum seal 
to be broken being that between the silica tube and the furnace endplate. This 
seal is made with vacuum wax. The lower endplate is similarly sealed to the 
silica tube, so that this can be evacuated. The viscosity determinations are carried 
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Figure 1. Figure 2. 


out in an atmosphere of hydrogen, which after purification by passing through 
chromous chloride and potassium hydroxide, is dried and admitted through the 
side tube P in the brass cover G, and allowed to escape through the outlet tube Q 
at the bottom of the tube. Before the hydrogen is admitted the furnace tube is 
evacuated through the tube O in the bottom endplate. 

In order to ensure that the phosphor bronze suspension wire be kept at a low 
constant temperature, and that the wax seals be kept cool, water cooling is provided 
at both ends of the silica tube. 

Measurements of amplitudes and times of oscillation are made by observing 
the graduated aluminium disc D through a telescope placed a few feet away from 
the apparatus. 
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§4. EXPERIMENTAL PROCEDURE 


A check determination was made of the viscosity of molten tin, the viscosity 
being calculated from the constants of the apparatus by the first of the methods 
described in § 2, that is without reference to observations on any other liquid. 
The value obtained was 0-01182 poise at 496:8°c. which agrees closely with the 
value 0:01187 obtained by Stott (1933). 

In order to be able to use the second method given in § 2, the apparatus was 
calibrated, using liquids of known viscosity and density. The results are given 
in the next section. In all experiments the depth of liquid in the crucible was 
maintained the same, the amount of liquid necessary to ensure this being weighed 
beforehand. 


The viscosity of molten zinc was then determined over a range of temperatures | 
between 430°c. and 480°c. The material used was high purity zinc, for which _ 
we are indebted to The National Smelting Co., Ltd. The analysis gave the | 


following quantities of impurities: Pb<0-001°% (probably 0-0005%), 


Cd 0-0001%, Cu 0-0001°%, Fe 0-0005°%, Caatrace. The metal was melted ina | 


separate furnace under a layer of wood charcoal, cooled and washed, first 
with 2%, acetic acid solution to remove surface oxide and then with boiling 
distilled water. It was then placed in the viscometer and the temperature raised 


to 480°c. After 24 hours at this temperature, the logarithmic decrement was 


determined and further readings taken 12 hours and 24 hours later. Constant 
readings having been obtained corresponding to 480°c., the temperature was 
lowered to 470°c., and after an interval of 24 hours at this temperature the 
logarithmic decrement was again determined. Measurements were made, 
following the same procedure at 460°, 450°, 440° and 430°c. ‘Trouble was at 
first experienced due to the vaporization of the zinc, but it was found that this 
could be eliminated by having a layer of powdered wood charcoal on the surface 
of the melt. Experiments on molten tin and molten zinc with and without such a 
layer of charcoal showed that the logarithmic decrement was not aftected by 
its presence. 


Suolon [8elehiey Oe ANS) 
(i) Calibration 


The constants in the equation 


2 3\ 1/2 
[E41 ]-anrsaepryssc (22) 
0 p 


were determined using five liquids: chloroform, water, carbon tetrachloride, 
benzene and tin, for which accurate values of 7 and p at various temperatures are 
available. ‘The relevant data are given in Table 1. The value of 5 given in the 
fourth column represents the observed value less the value of the logarithmic 
decrement obtained with the crucible empty at room temperature. This value 
was 0-003659 for the first four liquids and 0-001644 for the fifth. Each value 
given represents the mean of six separate observations. 

‘The values of 7 and p for the first four liquids given in the table are obtained 
from International Critical Tables. The values for tin are taken from Stott’s 
paper. Calculation of A, B and C from the above data by the method of least 
squares gives A = —()-0362, B=0-0438, C=0-0306. On using these values of 
A, B and C and the observed values of 5 to calculate the viscosities n of the 


The Determination of the Viscosity of Molten Metals FOL 


calibrating liquids from the calibration formula, the greatest discrepancy between 


calculated and actual values is found to be 0-0002 poise. The values of 8, p and T 
in the case of zinc (see below) are such that the determination of the viscosity 
using the calibration formula is an interpolation and not an extrapolation process. 


(u) The Viscosity of Zinc 


Using the above values of A, B and C the viscosity of zinc was determined 
at various temperatures. The results are given in Table 2. The logarithmic 
decrements 6 given in the table were obtained by subtracting from the observed 
values the logarithmic decrement with the crucible empty at 455°c. This value 
was 0001644. ‘The values of the density were taken from a paper by Pascal and 
Journiaux (1914). 


‘Lables i 
‘Temperature Viscosity 7 Density p Logarithmic sor 
Substance 5 : oscillation 
(aC) (poise) (gm/cm?*) decrement 6 T 
(sec.) 
Chloroform 23-00 0-005460 1:483347 0:007318 19-85 
Water 23-75 0-:009216 0:99738 0:007286 19-86 
Carbon 
tetrachloride 19-75 0:009743 1:594539 0-009678 19-90 
Benzene 19-75 0006495 0:879012 0:005692 19-87 
Tin 496°8 0-01187 6°785 0-023824 20-09 
Time of oscillation of empty crucible =19-81 seconds. 
Table 2 
‘Temperature Period Density Logarithmic Viscosity 7 
(C3) (seconds) (gm/cm?) decrement 6 (poise) 
480°8 20:26 6°833 0:03929 0:0378 
468-8 20-26 6847 0-03955 0-0383 
459-3 20-26 6-860 0-03981 0:0388 
448-7 20:28 6:876 0-04008 0-0393 
439°8 20-30 6-888 0:04027 0:0397 
429-6 20:31 6-903 0-04056 0-0403 


The relation between the viscosity 7 and the temperature is shown in the 


graph of Figure 2. 
Andrade (1934) has suggested a relationship between the viscosity of a liquid 
and temperature of the form 


Ue +/00 
HOW = Ae, 


where w is the specific volume, 6 the absolute temperature, and A and c are 
constants. On fitting this relation to the results for molten zinc it was found that 
with A=0-009533 and c=81-266 the calculated values of 7 agreed with the 
observed values to within 0:1%. 


(iit) Absolute Determinations 


The constants of the apparatus which enter into an absolute determination 
of the viscosity are the diameter of the crucible, the depth of the liquid in it, 
and the moment of inertia of the oscillating system. We have calculated the 
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viscosity of zinc using equation (16) and obtained the following values given in 
Table 3, the values of 5 and p being those of Table 2. 


Table 3 
Temperature (° c.) 480°8 468°8 459°3 448-7 439°8 429°6 
Viscosity (poise) 0:0380 0:0385 0:0390 0:0396 0:0400 0:0406 


The greatest difference between these values and those of Table 2 amounts to 3 in 
the last decimal place, i.e. about three-quarters of one per cent. 


(iv) Comparison with Previous Determinations 

The only previous work on the viscosity of zinc of which we are aware is that 
of Gering and Sauerwald (1935), who determined the coefficients of viscosity 
of a number of metals by observing the motion of small quantities of the metal 
along a capillary tube under the influence of gas pressure. They obtained values 
for zinc in the temperature range 450° to 700°c. ‘These values varied from 
0-03168 poise at the lower to 0-01865 poise at the higher temperature. The results 
of the present paper are higher than this, e.g. our value at 448-7°c. is 0:0393 
(Table 2), which may be compared with Gering and Sauerwald’s 0-03168 at 
450° c. 


§6. CONCLUSION 


The viscosity of pure zinc has been determined for various temperatures in 
the range 430° c to 480° c. by the oscillating cylinder method. It is intended to 
investigate the effect of adding small quantities of other metals such as aluminium 
to the zinc, and some preliminary measurements have already been carried out. 
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High Frequency Permeability of Ferromagnetic Materials 


By R. MILLERSHIP* anp F. V. WEBSTER 
Physics Department, University of Leeds 


Communicated by R. Whiddington; MS. received 20th March 1950 


) ABSTRACT. Methods are described for the measurement of resistive and inductive 


permeabilities over a range of frequencies from 150 to 10,000 Mc/s. The measurements of 
the resistive permeabilities were carried out by comparing the attenuation of a coaxial 
transmission line with the ferromagnetic as the inner conductor with that of a coaxial line 
having a non-ferromagnetic as inner conductor. The inductive permeabilities were obtained 
from the wavelength in a coaxial line having the ferromagnetic as the inner conductor. The 
results show that at all frequencies the resistive permeability is greater than the inductive, 
and both decrease to unity at very high frequencies. The effective permeability, which is 
assumed to be complex, is deduced as a function of resistive and inductive permeabilities and 
plotted against frequency for a variety of materials. 'The experimental values are compared 


~ with theoretical values derived from two formal treatments. 


Sele NE ROD UC LOIN, 
HE process of magnetization of ferromagnetics under the influence of 
ultra-high-frequency fields has until recently received very little attention, 
probably because of the difficulty of producing an oscillatory field of 
sufficiently high frequency, but, as a result of extensive advances in ultra-high- 
frequency technique during the war, this difficulty has been largely removed. 
Unfortunately most of the previous investigations have been confined to very 


) restricted frequency ranges; in fact, many have consisted of isolated determinations 
| of permeability at a single frequency. Furthermore, the materials used in the 


past have often been of ill-defined composition. As the low field properties of 
magnetic materials are so strongly structure sensitive, it is only by using the same 
material (and preferably the same specimen) over the whole frequency range of 
interest that the character of the variation of permeability with frequency can be 
determined with any certainty. 

A complete and critical account of the various methods which have been 
employed has been given by Allanson (1945). ‘The methods fall into two distinct 
groups: those in which the permeability is determined from a measurement of 
the resistive losses of a circuit containing the ferromagnetic material, and those 
which are dependent upon the effective reactance of a similar circuit. ‘The 
permeability deduced from measurements in the former category is usually 
denoted by jp, and that from measurements in the second category by p,. Of 
these two coefficients zp is always larger than p,, the difference between the two 
being covered formally by treating the permeability as a complex quantity. It 
isnecessary to determine both pp and py at the same frequency in order to calculate 
the complex permeability. This has not been generally realized, and, as a 
consequence, almost all the previous investigations have consisted of measurements 
of either jxp or jy, alone. 


* Now with the National Coal Board. 
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In order to provide conditions under which the field may be introduced into 
ne | 
ferromagnetic material. The parameters of a transmission line which contains 
ty are the attenuation constant, the Q factor, and the characteristic impedance, 


the specimen it is convenient to use a transmission-line system incorporating the 


whilst the wave velocity in the line gives the information necessary for the 
calculation of 4, Most of the measurements obtained by earlier workers were 
at frequencies of the order of 100 Mc/s., employing the open Lecher wire type 
of transmission line. Such a system would be unsuitable at higher frequencies 
as considerable power would be lost through radiation, and hand-capacity effects 
would also become significant. At these higher frequencies coaxial lines and 


waveguides have proved to be more satisfactory. One very accurate method — 


would involve the use of a waveguide constructed partly of the ferromagnetic 


material: the disadvantages of this are that the magnetic properties of the material _ 


may be affected by machining, and that the waveguide could only be used over a 
limited range of frequencies. 


One of the earlier theoretical treatments of magnetization at high frequencies _ 


| 


was that developed by Arkadiew (1926) in which the permeability was introduced — 
as a complex quantity. This theory contained much that is incompatible with — 


the modern theory of ferromagnetism, and a theory based on more modern 
concepts was developed by Becker (1938, 1939), who assumed that the decrease 
in permeability was caused by a reduction of the field in the ferromagnetic due 
to microscopic eddy currents induced by domain boundary movements. Kittel 
(1946) attributed the reduction in the permeability to a reduced effective field 
across the domain boundary caused by the small penetration of the applied field 
due to the skin effect. Earlier Landau and Lifshitz (1935), considering the 
change of magnetization as due to the change in direction of the magnetic vector 
in the crystalline anisotropic field, predicted the occurrence of a resonance effect 
when the frequency of the applied field is equal to the Larmor frequency for 
gyroscopic rotation of the electron spins about the crystalline field direction. 
In metallic ferromagnetics, whether the magnetization change is due to trans- 
lational boundary movements or to rotations of domain magnetization vectors, 
or both, the damping caused by microscopic eddy currents must be mainly 
responsible for the general decrease of permeability with frequency. These 
questions are further considered below (§ 4 (iv)). 


§2. METHOD 


(1) Measurement of the Resistive Permeability 


The method used for the measurement of pz has been described in detail 
in a previous paper (Hodsman, Eichholz and Millership 1949), and hence only 
a general outline will now be given. A coaxial-line system is employed with the 


specimen as the centre conductor. This arrangement has the advantage that — 


radiation can be propagated down the line in the fundamental mode over an 
infinite range of frequencies and so the introduction of the oscillatory field presents 
little difficulty. As the specimen is in the form of a wire, accurate machining to 
conform with the geometry of the system is eliminated. Furthermore, as the 
outer conductor is earthed, radiation is prevented from entering and leaving the 
system, which is consequently immune from external interference and instability 
due to hand-capacity effects. 
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The permeability is determined from the attenuation constant of the system, 
which is given by (see, for example, Jackson 1945, p. 50) 


9-95 x 10-8 F271 fu \N2 1/4 \¥? 
— log (b/a) _ Re a (2) | ab/mette eee: (1) 


where a, 4,, and o,, are the radius, permeability, and conductivity respectively 
of the inner conductor, 6, 44, and o, are those of the outer conductor, and f is the 
frequency of oscillation in megacycles per second. If a line of length J is used 
and the attenuation is measured with the inner conductor firstly of ferromagnetic 
wire and secondly of non-ferromagnetic wire, then the difference in attenuation 
) 1s given by 


Ben tila, \i= a | (2)" . (-)" | dbomeeoes (2) 


where jz, and 9, refer to the ferromagnetic specimen and o, to the non-ferromagnetic. 
' In this way the effect of the outer conductor and of any other residual attenuation 
is eliminated. Measurements on the standing wave pattern in the line are trans- 
formed into values of input impedance for different terminations and plotted as 
circle diagrams of impedance, from which the attenuation constants are derived. 


(11) Measurement of the Inductive Permeability 


It is apparent from the results of previous investigators that the dispersion 
in pz, takes place at frequencies considerably below those at which the change 
iN fp is most rapid. Preliminary investigations indicated that, at 1,500 Mc/s., 
pz, Was equal to unity for all the materials employed. Apparatus for the measure- 
ment of 4; was therefore designed to cover the range from 150 to 1,500 Mc/s. 
only. 

In order to measure yz, under the same conditions as pp, a section of coaxial 
© line was used in preference to the more easily constructed Lecher line, which 
4 ensured that the field configurations were precisely the same in both cases. ‘The 
+ method employed is similar to that described by Glathart (1939), being based 
i on the measurement of the wavelength in a coaxial line having the specimen as 
» its centre conductor. Writing d, and d for the free space and measured half- 
- wavelengths respectively the permeability is expressed as (Glathart 1939, p. 834) 


. a b\2 (d,? — d?) < 
py=l 256 (<) (Iog ‘) i dee x 10 Es attecniielrs (3) 


i where p is the resistivity of the inner conductor. Equation (3) may be solved 
| for Ad giving 


A. £ 3-544 x 10%a log (b/a) \ ‘| 
Ad=dy—d=d,| 1 {ape Se Totton BV) , 


By suitable choice of the dimensions of the system, the half-wavelength 
' shortening d)—d may be made sufficiently large to allow an accurate practical 
+ determination of p;, to be made; this is so when a0, and ba. Having 
- regard to mechanical stability, accurate centring of the inner conductor and 
» convenient operating conditions, a fairly large diameter outer tube containing a 
' very fine inner conductor wire was designed. A schematic diagram of the 
| apparatus is shown in Figure 1. 
| PROC. PHYS. SOC. LXIII, 10—B 3H 
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seen momar 


The wavelength is determined by altering the length af the line by means att 
a short-circuiting piston until two successive resonance peaks have been located, 
The distance moved by the piston between the resonance positions is then half 
the wavelength in the line. Resonance is observed by extracting a fraction of 
the power in the line by means of a small loop introduced through a slot at the 
input end of the outer conductor. The high-frequency current induced in the 
loop is rectified by a crystal and the output displayed directly onameter. Figure2 
shows some typical resonance curves. 


Resonator 


Oscillator 
Figure 1. Layout of apparatus for measurement of py. 
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Figure 2. Typical resonance curves. 


For the calculation of the free space half-wavelength equation (4) may be 


written 
Ng erate sae 
is =i (Gaye) |: ae ©) 


where Ad, is the difference between the free space half-wavelength and 
the measured half-wavelength for a non-magnetic inner conductor and | 
k =3-544-x 10%a log (b/a). As Ady is always small, dy in the above equation may | 
be replaced by the measured half-wavelength. 

The free space half-wavelength is then obtained as the sum of Ad, and the 
measured half-wavelength for a non-magnetic inner conductor, i.e. 


dy=d |2 s (arc) | pn Es ee (6) 
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$3. RESULTS 

The following materials have been investigated, all in wire form: Swedish 
charcoal iron, white annealed and bright; 50/55 carbon steel, white annealed 
and bright; mumetal 6, hard drawn. 

The percentage compositions and treatments of the specimens were as follows : 

Iron: C, 0-03; Mn, 0:05; S, 0-002; Si, P, Cr, Ni, Co negligible. 

Steel: C, 0:50-0:55; Mn, 0-5-0-6; S,0-03; Si,0-1; P,0-02; Ni,0-2; Crnil. 

maumetal> Ni, 77; Fe, 153 Cu, 5; Cr, 2; Mn, 0°8; Si, 0-1. 

The iron and steel specimens were hard drawn with 90° cold work. The 
white annealed specimens were annealed with a wet surface at a temperature of 
700° c. in an inert atmosphere of town gas (CO, CO,, N,) and the bright specimens 
were left unannealed. 
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Figure 3. Relative resistive permeabilities. 
Curve 1: Steel bright. Curve 2: Iron bright. Curve 3: Iron annealed. Curve 4: Mumetal. 
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Figure 4. Relative inductive permeabilities. Figure 5. Relative inductive 
Curve 1: Steel bright. permeabilities. 
Curve 2: Steel annealed. Curve 1: Iron bright. 


Curve 2: Mumetal. 
Curve 3: Iron annealed. 
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In order to provide a firmer basis for comparison the results are shown as 
) relative values, i.e. p/w and pyz/u, in Figures 3, 4 and 5. ‘The values of the 
initial permeability , are: iron annealed, 180; iron bright, 102; steel annealed, 
103; steel bright, 65; mumetal, 440. 
| At any given frequency jp is greater than ju, for all the materials employed. 
_ Both pp and py, appear to decrease to unity, although the frequency range 
3 H-2 
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employed was not sufficiently great to obtain such a value for py. No values . 
below unity have been obtained, such as those reported by Birks (1948) for finely 
divided ferromagnetic compounds in paraffin wax. Owing to the lack of suitable 
oscillators it has not, as yet, been possible to obtain a sufficient number of readings ; 
to cover adequately the range where j1, approaches unity. 


§4. DISCUSSION OF RESULTS 
(i) Introduction 
An estimate of the strength of the high-frequency field that is applied to the 
specimen may be obtained from the equation for the power flowing down the » 
line, 


P=n (es) H,2log(b/a),, ee (7) 


e 


which is derived by integrating the time average Poynting vector over the cross | 
section of the line. Assuming all the power produced by the oscillator is dissipated | 
in the line an upper limit of 0-1 oersted is obtained for the field strength. The 
experimental values for the permeability therefore refer to magnetization in the 
reversible pcrtion at the commencement of the magnetization cycle. 

Any theory of the effects must explain the general decrease in the permeability, 
and also the difference which exists in the values obtained by resistive and inductive » 
experimental methods. 


(ii) Effect of the Surface Condition 

It is apparent that the permeability for a high-frequency field is a property 
of the surface layer of the material and consequently will depend to some extent 
on the nature of the surface. It has in fact been suggested by Wien (1931), 
Michels (1931), and Procopiu and d’Albon (1937) that the general decrease in 
permeability with increasing frequency of the field is due to the existence of a — 
thin surface layer (35-70 uu) of non-ferromagnetic material. Since the thinnest — 
truly ferromagnetic films obtained by Sorenson (1924) and Tyndall (1924) | 
are of the order of 20 wp, the existence of such a non-ferromagnetic layer appears | 
to be unlikely. It is probable, however, that the domain structure in the surface | 
layer will be different from that in the main body of the material : hence, if boundary ~ 
movements are responsible for the magnetization at high frequencies, the magnetic - 
properties will vary considerably for different skin depths and therefore for different — 
frequencies, 

Unfortunately, the permeability and conductivity always appear together as | 
a ratio, i.e. y/o, in the equations for the attenuation and wavelength in a coaxial | 
line. It is known that a change in conductivity takes place for non-ferromagnetics | 
when the surface is etched or scratched and it is expected that a similar change | 
takes placeforferromagnetics. Ifthesurface is roughened the increased conducting © 
path gives rise to a decreased conductance and permeance. These quantities are | 
both dependent on the length of the conducting path and are reduced equally. 
Due to the reduction in the effective values of conductance and permeance the 
skin depth increases and brings about a corresponding increase in these quantities, 
the ratio y/o in the equations always remaining the same for different surface 
conditions. ‘Thus it is very difficult to relate the permeability alone with surface 


condition, although it is intended to carry out experiments in this field in the 
future. 
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(iit) Difference between frp and uy, 


The most satisfactory approach to this problem was proposed by Arkadiew 
(1932) who treated the permeability as a complex quantity, i.e. of the form 
#@=fy—tMy,. Considering a wire of ferromagnetic material, wp may be regarded 
as that value of « which satisfies the theoretical equation for the resistance. 
Similarly yu; emerges from the equation for the inductance. 

The high-frequency resistance of an isolated conductor of radius 7) per unit 
length is given by Jackson (1945, p. 38) as 


1 (afp\? 
= Dut (#) fo eo St (8) 
which gives (in C.G.s. units) 
R 7 (27wpx\"? 
Rer2s ( ; Sey ea ee (9) 
where R, is the D.c. resistance of the wire and p is the resistivity. Similarly 
wl 7o (27wpy,\1? 
Rae a P cee Sate aC ee Pecan (10) 
Maxwell’s equations for the specimen are 
OH _tE ) 
ara cp | 
ae ee) g gets (11) 
Ok _ ov 
or S Cc i 


where 7 is the radial axis perpendicular to the direction of the wire. 
Eliminating H, 


gd OS a Ne mg 
Sn ed ee ee 2 
aye = 1 ep E (12) 

Substituting a complex permeability, 4 =p, —7pg, gives 
4 eae bee 
E=Cexp [| a Getta) r|, a ee (13) 
OE 47w : te 

~— =\(—= ‘gage Es See ae se 14 
and ay { ep (He + i)| E (14) 


The effective resistance R and inductance L are obtained from (Abraham and 
Becker 1932, p. 199) 


Rtiwl 15 ‘at ( E ) ye (15) 


ie <i cp \dElor 
= 2 {ae (Hebin)} nee (16) 
= F(A) ant tate oad LG + 18)! Doo) 
and therefore Me (Met oe?) + oe aes er Wabi veer x (18) 


a) 


ain (14? + Mg — Fe 
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from which Hy =(HRer)? | 
| 


H2=3(MR-PL) J 
These equations enable us to calculate 4, and p, from the experimental values 
of wp and py. P 


(iv) Theoretical Treatments : : 
(a) Becker theory. 


Becker (1938) considered the damping effect of microscopic eddy currents, | 


associated with the movement of the domain walls, to be responsible for the 
general decrease in permeability with increasing frequency. As the domain wall 
moves, eddy currents, which are proportional to the velocity of movement, are 
induced in the surrounding medium. ‘The eddy current field will be opposed 
to the change in magnetization and may thus be interpreted as a braking field 
which reduces the effective field acting on the wall. ‘The displacements will 
therefore be reduced, resulting in a reduction of the permeability. This treatment 
leads to an expression of the form 
fo— 1 
pw—-l= idol)... (20) 
where jy is the D.c. initial permeability and w, is a ‘critical frequency’ which, 
for a cubic domain of side / embedded in a ferromagnetic medium of susceptibility 
Ko, 1S expressed as 
302 
= 4moxl? Cin SR De ee ea Dee 


The high-frequency field is considered to be constant over the whole volume of 
the domain; owing to the small skin depth, however, the applied field may not 
even penetrate to the extent of the individual domains in the surface layer for 
centimetre waves. As the domain model assumed is over-simplified the equation 
for w, is of doubtful validity, although equation (20) will still hold. Such an 
expression will result from any consideration of a damping force which is 
proportional to the velocity of movement of the boundary; and it is reasonable to 
suppose that these damping forces arise from the action of microscopic eddy 
currents. 


Separating the real and imaginary parts of equation (20) we get 


Wy 


ta Ho 1 ] 
1+ (w/a,)? | 


_ (Ho=Wlwlor) | 
mo T+(wjnPe | 


the equations for the components of the complex permeability. 


(b) Kittel theory. 


i= 
aco Ce 


A theory developed by Kittel (1946) relies on a more thorough consideration — 
of the skin effect in a film one domain thick. The effective field acting on the — 


domain boundary is the mean value of the applied field taken over the entire 
wall which moves as a rigid whole under its influence. The reduced permeability 
is thus accounted for in terms of the small effective field. 

The theoretical curves are similar to those derived by Becker for frequencies 
beyond 3,000 Mc/s. but at lower frequencies they are very much steeper, and a 
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} resonance peak appears on the initial part of the wp curve: of this there is no 
“| suggestion in the experimental curves. It seems therefore that the reduced 
penetration of the field due to the skin effect cannot be the sole cause of the 
ferromagnetic dispersion. 


( (ce) Damped resonance considerations. 


The results of experiments performed by Griffiths (1946) and Birks (1950), 
» with a static field imposed on the specimen in addition to the transverse high- 
+ frequency field, indicate the existence of damped resonance phenomena. 
Resonance due to gyroscopic rotation of the electron spins about the field direction 
may occur when the frequency of the oscillating field is equal to the natural 
} frequency of gyration given by the Larmor precession frequency, wy =eH/47me, 
where H is the strength of the field along the axis of rotation. Withasuperimposed 
static field the electron spins may be aligned throughout the material and a sharp 
» resonance condition is obtained. A fuller treatment has been given by Kittel 
‘| (1947, 1948) who showed that the resonance frequency depends on the shape of 
) the specimen, and fora plane surface is given by w,) =)(BH)"? (rather than w) =yH), 
ay where y = e/47mc =2-80 Mc/s/oersted. 

\ Landau and Lifshitz (1935) were the first to consider resonance in the absence 
of an external biasing field, and deduced that it would occur about the anisotropy 
» field direction. In an unmagnetized polycrystalline material the internal field 
will vary in direction throughout the volume and if strains are present there will 
also be a variation in magnitude; hence each domain will resonate at a different 
frequency and a broad flat resonance will be obtained instead of a sharp one. 

In the experiments under discussion no static field was applied to the specimen 
and so the experimental results may be interpreted in terms of a damped resonance 
occurring about the crystal anisotropy field direction, the damping being caused by 
microscopic eddy currents. 

The resonance dispersion equation for an oscillatory field transverse to the 
internal field is given by (Frenkel 1945) 


poole” Co" 


Byo—-1 wy? — w? 4+ Ziwe,’ 
where w, is the ‘resonance’ frequency and w, is the ‘damping’ frequency. 
Separating the real and imaginary parts of equation (23) we get 

pio Wo" (Wo — w”) 
Mo—1 (wo? — w?)? + 4w2w,? 


[Hy 20970 


for the components of the complex permeability. 


(d) Summary of theories. 

The three methods of attack presented above may be summarized as follows: 
(1) Becker theory. ‘This theory is based entirely on damping produced by 
eddy currents on domain boundary movements. ‘The formal treatment gives 
for ,. the equation 


__ Mon! 
EaeesiceiGaian)’ ce eunreld (20) 
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which contains one parameter that can be adjusted to give as close a fit to the 
experimental curves as possible. 


(2) Kittel theory. Damping is not taken into account in this theory, the 


dispersion being wholly attributed to a reduction of the effective field acting on 
the domain wall because of the very small skin depth. 


(3) Damped resonance. Rotations only are considered here and the formal | 


treatment gives for ju 


a (23) 


Mo—1 a — w+ 2iww,’ 


an equation with two parameters, both of which are adjusted to give the best fit 
to the experimental curves. It seems probable that these effects are present, 
not alone but all together in varying degrees, and, therefore, it is not to be expected 


that a perfect fit to the experimental curves will be obtained by considering any | 


one effect. The aim is to obtain a general fit which may give some indication 


of the type of mechanism responsible for the magnetization at high frequencies, — 
e.g. rotational movements of the spins or translational movements of the domain | 


walls. 


(v) Comparison of the Experimental Results with Theories 
(a) Becker theory. 


A comparison of the experimental results for the complex permeability with 
values of 4, and py, calculated from equation (22) are given in Figures 6, 7 and 8. 

The theoretical curves are fitted to the experimental results by making w, 
(in equation (22)) approximately equal to the frequency at which the experimentat 
value of jz. reaches a maximum. From the values of w, the sizes of the domains 
are calculated using equation (21). These are: iron annealed/=1-7 x 10-*cm. ; steel 
annealed /= 1-4 x 10-4 cm.; mumetal /=1-6 x 10-4 cm. These values are of the 


right order of magnitude although little reliance can be placed on them because of — 


the over-simplification of the domain model. It has been pointed out by Landau 
and Lifshitz (1935) and Néel (1944) that the domains in the body of the material 
will most likely be in the form of plane parallel sheets, whereas the surface layer 
will consist of small domains of triangular cross section which provide ‘ closures 
of flux’ for the main internal domains and thus prevent the appearance of free 
poles at the surface. 

Although the general forms of the experimental and theoretical curves are 
similar, the frequency at which p, is equal to jy is in all cases lower than that at 
which jz, is a maximum, whereas the theory shows that the ,., curve should pass 


through the maximum of the p, curve. Also, the real part of the permeability — 
decreases more rapidly than the theory indicates; this was observed by Johnson | 


and Rado (1949) who experimented at frequencies of 200 and 975 Mc/s. with a 
superimposed static field applied to the specimen. The theory gives a much 
larger and sharper maximum for y. than the experimental curve suggests. 


(b) Damped resonance. 


‘Theoretical values of the complex permeability are derived from equation (24) 


and are plotted against frequency together with experimental values in Figures 9, 
10 and 11. 


a 
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The curves are fitted to the experimental values by choosing w, and w, to 
produce a uy curve which has a maximum approximating to the experimental 
maximum. ‘The condition for a maximum in py is 


SRE S20, Nosh ee ea (25) 
10 Lon 
[ L,-| Experimental 
Wel —— Experimenta Ral —-—— Theoretical 
08 —-—— Theoretical OS a Poo! (@,= 500 Mt/s.) 


\ (w= 350 Mc/s.) 


0-6 


=< 


04 0-4 
Ree Ss. +k he 02 
iy l ae | ==] 
CF 500 1000 1500 2000 °F 500 1000 1500 7000 
Frequency (Mc/s.) Frequency (Mc/s) 
Figure 6. Comparison of experimental Figure 7. Comparison of experimental 


results with the Becker theory— 
lron annealed. 


results with the Becker theory— 
Steel annealed. 
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Figure 8. Comparison of experimental results Figure 9. Comparison of experimental results 
with the Becker theory——-Mumetal. with the damped resonance theory—Iron 
annealed. 
10 
a= — A plete! Experimental 
08 tym ---—-— Theoretica ----— Theoretical 
\. (,=1950Mo/s. and w,=4000 Mc/s) (= 1500 NMo/s. and w,=4000 Majs) 
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i Frequency (Me/s.) lrcauency (Me/s) 
; ; ; ‘ 
Figure 10. Comparison of experimental results Figure 11. Comparison of experimental results 
| with the damped resonance theory—Steel with the damped resonance thcory— 
i annealed. Mumetal. 


Due to a lack of suitable oscillators covering the range of frequencies between 
» 500 and 1,000 Mc/s. the exact location of the maximum on the experimental curve 
for jz. cannot be obtained without further results for y, in this range. 
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In some respects the type of curve obtained from considerations of a damped 
resonance expression approximates more closely to the experimental results than 
that drawn from an expression of the Becker type. The , and ps, theoretical 
curves intersect at a frequency slightly lower than that at which py, is a maximum, 
which is in better agreement with experiment. Also, the , curve is nearer to 
the experimental results than the Becker pz, curve. 


There are certain points, however, in which the experimental and theoretical — 


curves do not agree. Firstly, the theory indicates that 4, becomes unity at a 
frequency equal to w, and becomes less than unity at frequencies greater than wo. 
Secondly, the theoretical curve (like the Becker curve) has a much larger and 
sharper maximum than the experimental points suggest. . 

A comparison of the resonance frequency with that which would be obtained 
for a resonance of the spins about the anisotropy field may be made by substituting 
the value of the field in the Larmor expression. The anisotropy field is given 
(Birks 1950) by H=2K/I,, where K is the anisotropy constant and J, is the 
saturation magnetization. For iron, this is approximately 500 oersteds, which 
would give a Larmor frequency of 1,400 Mc/s. The resonance frequency for 
iron obtained from the experimental curves is therefore of the same order of 
magnitude as the Larmor frequency for the electron spins about the anisotropy 
field. ‘The Larmor frequencies for steel and mumetal are approximately 1,400 
and 120 Mc/s. respectively, which, especially in the case of mumetal, is not in 
very good agreement with the value of wy substituted in the equation for». Itis 
unlikely that an entirely different mechanism is operative for mumetal, and seems 
to indicate that more than one process is involved, e.g. a combination of boundary 
movements and spin rotations, both being damped by microscopic eddy currents. 


y5., CON CLUSION 

Experimental values of the permeability have been obtained from both 
resistive and inductive measurements on a circuit containing the ferromagnetic. 
The values of np and pz; both decrease to unity as the frequency of the applied 
field is increased and jp is always greater than ju,, the difference between them 
being covered formally by supposing the permeability to be complex. 

The experimental values of the components of the complex permeability are 
compared with values obtained from two theoretical treatments, one based on 
boundary movements and the other on rotations; a fairly good fit is obtained in 
each case. ‘The resonance frequency for iron calculated from damped resonance 
considerations approximates to the Larmor precession frequency of the electron 
spins about the anisotropy field, but the divergence for mumetal suggests the 
existence of more than one mechanism. It seems probable that the magnetic 
dispersion is due to eddy current damping of translational and rotational move- 
ments to approximately the same degree, the situation being further complicated 
by the possibility that skin depth effects of the type suggested by Kittel are 
responsible in some measure for the dispersion. 
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The Production of Pulsed Magnetic Fields 
using Condenser Energy Storage 


By Kees...W. CHAMPION 


Electron Physics Department, University of Birmingham 


Communicated by F. Sayers; MS. received 21st March 1950 


ABSTRACT. Aneconomical way of producing a pulsed magnetic field is to use condensers 
as the energy source and an air-cored magnet coil which, with the storage condensers, forms 
an oscillatory circuit. By means of a switch the coil is connected to the charged condensers 
and then disconnected precisely at the end of one cycle. In this way all the stored energy is 
used to produce the field, while most of it is recovered at the end of the pulse, the only 
appreciable losses being due to the resistance of the coil and switch. 

The design of coils for use in the way described above is considered in detail. Expres- 
sions and curves are presented which enable the field produced by a specified coil and 
condenser combination to be readily calculated. Conversely, the data can be used to design 
the most satisfactory coil for any practical application. An example is given of the design 
and construction of an air-cooled coil. 


Sie UNC RODMIC TO IN 
O investigate certain fundamental properties of low pressure arcs it became 
| necessary to study the effect on them of strong magnetic fields. As the 
arcs in question were to be pulsed and the field required was inconveniently 
large for a conventional iron electromagnet, it was decided to pulse the magnetic 
field. The most satisfactory sources of the current pulse for such fields are 
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suitable generators or condensers. The design of coils for use with generators — 
has been considered very fully by Cockcroft (1928). 

However, very often condensers are a more satisfactory source of energy — 
when such considerations as maintenance are taken into account, and, in the 
following, the design of suitable coils for use with a condenser energy source is 
considered. It is hoped that these data will be of use when it is necessary to 
design efficient coils for producing strong pulsed magnetic fields, for example, 
those for use with cloud chambers. 


§2. THEORY OF PULSED MAGNETIC FIELD 
If it is not necessary that the magnetic field pulse be square, a convenient 
and economical way of producing it is as follows. With the circuit as in Figure 1 _ 
the switch is closed when a pulse is desired and opened again precisely at the end — 
of one cycle. By this means all the stored energy is used to produce the magnetic — 
field, while most of it is returned to the condenser at the end of the cycle, the | 
only appreciable losses being due to the resistance of the coil and switch. In > 


practice it is usually found that an electronic switch is much more satisfactory — 
than a mechanical switch. | 


Inductance L eee Fs 


Se 


WLLL 


All lengths in cm. 


Charging Circuit 


Current 
v incoil 


Voltage. 
across Coil 


Figure 1. If R is not negligible the phase of the Figure 2. 
current is shifted and the waves are damped. 
If it were possible to convert all the energy stored in the condenser into a 
uniform magnetic field Hy within the cylinder enclosed by the inner diameter 
of the coil, then the following relation would hold 


He 


ina ma,72l = ACV? ' 10 
k 20C\1l2 V 
1.e H,= (=) a, aL oroePane (1) 


where /) is in gauss, Cin wr. and V in volts and it is assumed that the permeability 
is unity. _‘I’he symbols representing the coil dimensions are defined in Figure 2. 
Thus the actual magnetic field H at the centre of a coil can be expressed by 
20C\12 V 
H=(——) —kK 
( i ) a K one (2) 
where K is a positive quantity less than one. It is convenient to split K into 
the product of two quantities S and J where S depends only on the shape of the 


i 
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coil and arises from the fact that the field is not uniform and not all contained 
in the cylinder enclosed by the inner diameter of the coil. J, on the other hand, 
is a factor to correct for the reduction in the first current peak (and, hence, the 
maximum magnetic field) due to the finite resistance of the coil. The values of 
S and J each lie between zero and one. 
He (AF)"s cy ec ae Pesce ove 3) 
i a, 
The value of the magnetic field at the centre of a coil may be expressed in the 
form 
2rnI Ay + (/? + a,?)¥2 
10(a,— 4) © is a, + (P+ ay?) 
where 7 is the number of turns in the coil (assumed to be uniformly distributed) 


H= 
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Figure 3, L=n’?a,A or n=(L/a,A)!” L in pH., ad, in cm. 


and J is the current in amperes. Now n=(L/a,A)"?, as indicated by the curves 
for inductance* (Figure 3), and the maximum value of J is 

Ve-RT/8L 
(R? + L?w?)"2? 


1 RR 2 
where w= Lc = 412 and T= mt é 


l= 


* These curves were calculated from data in a paper by Grover 1922. 
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Figure 4. 
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If we substitute for m and J in (4) we obtain 
2r LAV JCv2 a, + (1? + a,?)W? 


lak CREE A RUD ALP eM CEE 
From comparison of this expression and (3) 
1/2 2 2\ 1/2 
a ees a sclenaes 
iL 1 
Seren ea ae ene ee (6) 
10o<52* (2 -1)," 1+ (= +1) 
a ay 


This is the required shape factor since it is a function of //a, and a,/a, only. 
Curves of S are plotted in Figure 4. As would be expected, the value of S 


tends to unity for long, thin coils. 
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Figure 5. 
R p10° R resistance of coil in ohms., p specific resistance (ohm. cm.), 
L a,2As¢ L inductance of coil in henries, a, internal coil radius (cm.), As space factor. 


Now R/L=10%/a,?A,6 where p is the specific resistance of the material of 
the coil (in ohm.cm.), A, the space factor* and ¢ a quantity which is a function 
only of the shape of the coil. 
2 ON SITAR 7) 
=— —A{———_}. 
$ 7 ag ay Gres 


* The proportion of the coil winding volume occupied by conducting material is called the 


space factor. 


800 K. S. W. Champion 


Values of 4(f) are plotted in Figure 5 and values of J as a function of RT/L are 
set out in the Table. For values of RT/L up to 1-8, J may be represented by 
(1— RT/8L) with an accuracy closer than 0:2°%. 


Calculation of K for a Coil 


Thus, to calculate the magnetic field produced at the centre of a coil, operated | 
in the manner previously described, it is necessary to calculate K. This can 
be done by making use of the curves in Figures 4 and 5 and the values in the 
Table. The value of K is completely determined by J, a, a, A,, p and T. 

As an example, consider a coil with ]=5 cm., a, =2-5 cm., a, =4 cm., A, =0-7, 
p=1-:724 x 10-* ohm.cm. and T=1 millisec. In this case l/a,=2 and a,/a,=1-6 
and thus ¢=3-15 and R/L=125-1. Now RT/L=0-1251 and hence J =0-984; 
since S =0-803 the value of K is 0-790. 


‘Table 
Difference per Difference per 
RE = 0-1 of RTL Be z! 0-1 of RTL 
0 1 Ded 0-693 0:0113 
0-01 0-999 3 0-638 0-011 
0:05 0-994 3355) 0-586 0-0104 
0-1 0-987 4 0-537 0-0098 
0:2 0:975 Per 1 of RT/L 
0-3 0-962 5 0-451 0-086 
0-4 0-950 6 0-379 0-072 
0:5 0-937 0:0125 ey 0-319 0-060 
0-6 0-925 8 0-270 0-049 
0:8 0-900 10 0-195 0-038 
1:0 0-875 12 0-143 0-026 
ily 0-850 15 0-0923 0:0169 
1:4 0-825 20 0-0463 0-0092 
1-6 0-800 25 0-0239 0-0045 
1:8 0:775 30 0-0125 0-0023 
2:0 Oe75i 0-012 40 0:0035 0-0009 
2-2 0-727 0-012 50 0-00099 


ioe) 0 


If the coil is used with a condenser bank containing 1,000 joules a field of | 
19,980 (approximately 20,000) gauss will be produced at its centre. Note that, 
with the above specifications, the ratio R/L is determined. Also, since 

1 R2\-12 
LC 41? 
the product LC is determined. Hence, when any one of R, L, C is specified the 
values of the other two are fixed. 


T=2n( 


Limitations of the Data 
Before considering the design of coils it is desirable to point out the limitations 
of the application of the data contained in this paper. In the curves for inductance | 
(Figure 3) it has been assumed that the inductance is independent of the space 
factor, which is true when the latter is near unity. According to Grover (1922, 
p. 462), “only when the wires are relatively far apart will the correction (to be 


| The quantity ¢ is identical with Cockcroft’s ¢, (1928). However, additional values have been 
calculated for the curves in Figure 5. 
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} added to the inductance) exceed a few tenths of 1%’. As these inductance 
+) curves have been used in the calculation of the ¢ and S curves any deviation 
» from the former will be reflected in the latter. Secondly, the data only apply 


accurately to relatively low frequencies, as at radio frequencies the calculated 


rf inductance is reduced slightly and the resistance increased by the skin effect. 
) Also, at these frequencies the coil capacity becomes appreciable and radiation 


losses may no longer be negligible. 


$35 EXAMPLE OF THE DESIGN OF A COIL 
As a source of energy a condenser of capacity 7 uF. and maximum working 


_ potential 25 kv. is available. Consider the design of a coil to produce a maximum 


magnetic field with this source, given that its length is to be 20 cm. and internal 
radius2-9cm. ‘The pulse length is to be 10 milliseconds and, from considerations 
of the insulation requirements, a space factor* of 0-7 seems possible. Assume 
that the mean temperature of the coil is 20°c., then the value of p is the specific 
resistance of copper at this temperature. Values of K calculated from these 
data are plotted as a function of a,/a, in Figure 7. K has a maximum value of 
0-759 when a,/a, =1°5. 
Heating of Coil 

At this point it is necessary to consider the energy dissipated in the coil and 

the resultant temperature rise per pulse. ‘The temperature rise per pulse 


T Rj2 
t= | RidT 


0 vod 


where z is the current at any instant, v the volume of the coil windings and d and o 
respectively the density and specific heat of the material of the coil. Assume 
that R is constant during a pulse and take $/,? as the mean value of 2? where J, 
is the amplitude of the current vector at $7. ‘Thus, with sufficient accuracy 


me 
~ Quod 


where J, =VJ,(C/L)"? with J; =1—RT/4L for small values of RT/L. Values 
of 5¢ calculated from this expression are plotted in Figure 7. 

Now it is desired to operate this coil with one pulse per second for considerable 
periods. With a,/a,=1-5 and neglecting thermal conduction the rise in tempera- 
ture per minute of operation would be 32-7°c. Even with mica insulation and 
a Tufnol former it is not advisable to have an increase in temperature greater 
than 150°c. and for this reason, a larger coil is required. For example, with 
@,/a, =2-4 the temperature rise is only 4-1°c. per minute. Such a coil, even with 
no loss of heat, could be operated for at least 30 minutes and, with efficient cooling 
of the outer surface, this time could be extended indefinitely. 

Choosing a coil of this size the appropriate wire gauge is now to be determined. 
The resistance per unit length 


Ra}??? 
7 L(y + ay) 


bt 


= 1-074 x 10-4ohm/cm. 


* A higher space factor could be achieved if wire of rectangular section were used. 
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Hence 17 s.w.c. is the nearest gauge as it has a resistance of 1-068 x 10-4 ohm/cm, 
If wire of the calculated resistance were used 3,520 turns would be required. 

A coil based on this design has been constructed and is illustrated in Figure 6. | 
There are only 3,410 turns on the coil as its length is slightly less than 


Figure 6. Coil in its mounting. The coil can be supported with its axis either vertical 
or horizontal by screwing the base plate to the appropriate part of the frame. 
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Note. In this Figure T=10 msec. should read T=10-7 msec. 
2) cm, and 17s.w.c. is very slightly thicker than required. With this coil 
H =21,200 gauss with C=7 ur. and V =25,000 volts.* 


* Since the coil was constructed the condenser bank has been increased to 8 pP., giving 
H=22,700 gauss and T=10:7 msec. 
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Mechanical Stresses in Coil 


With intense magnetic fields the stresses in a coil may become very large. 
Cockcroft (1928) has given formulae and data which enable these stresses to be 
calculated. ‘The main resultant force is perpendicular to the axis and outwards 
(compensated by binding the coil) but with very intense fields (say, of the order 
of 500 kilogauss) the forces on the inner parts of the coil may be so great as to 
cause the copper to flow towards the outer surface. 

The maximum values of the stresses occur at the mid-point of the external 
surface and, with the coil designed above, are 6:9 kg/cm? in a direction perpen- 


) dicular to the axis and 4-3 kg/cm? parallel to it. ‘The total bursting force would 


be approximately 4,800 kg. ‘The copper windings should stand this force, but, 
nevertheless, the coil is reinforced with eight }-inch thick radial Tufnol strips, 
securely held by the end cheeks, pressed against the outer surface of the coil 
and with a brass band round their outer edges. The bolts securing the brass 
band are insulated from one end to minimize eddy currents. ‘The ‘Tufnol strips 
are set radially to facilitate efficient cooling of the surface of the coil. 


84. Kk AS A FUNCTION OF ITS PARAMETERS 


K may be regarded as an efficiency constant for the operation of a coil in the 
way specified at the beginning of this article. However, it must be pointed out 
that the field produced also depends on the energy stored and the volume enclosed 
by the inner diameter of the coil. 


0-7 = 
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0-6 = A a,=lcem.  y5=07 
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a,/a, 
Figure 8. 


Figure 8 gives examples of the way in which K (as a function of a,/a,) depends 
onJl/a,. The value of a,/a, for maximum K does not change much but, as would 
be expected, K is much smaller for shorter coils. In Figure 9 can be seen the 
very great change in K produced by changes in a,. This arises from the fact 
that, with smaller a,, the resistance of the coil becomes much more appreciable 
and, not only is K much smaller, but the value of a,/a, at which the maximum of K 
occurs is larger. Figure 10 illustrates the way in which K depends on T. 

31-2 
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Now, the only way in which K depends on parameters other than those of — 

shape (i.e. //a, and a,/a;) is in the expression for RT/L (of which J is a function) ; : 

RE pled | 

EL apPrd- 

The effect due to a, varies inversely as its square, whereas that due to 7 is — 
proportional] to its first power. From the expression for RT/L we see that 
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any proportional change in p or 1/A, has the same effect on the value of K as the 
same proportional change in 7. 

The curves in Figures 8, 9, 10 and the above considerations should be of use 
when a coil is to be designed, making it possible to assess rapidly the influence 
of the various parameters on the magnetic field produced. 


Production of Pulsed Magnetic Fields 805 


§5. SPATIAL VARIATIONS OF THE MAGNETIC FIELD 


In many investigations one other property is important, viz. the degree of 
» uniformity of the magnetic field. As the requirements of uniformity depend 
+ very much on the particular investigation, no general criterion will be presented, 
but merely data which make it possible to calculate the degree of uniformity of 
1 the field produced by a particular coil. As would be expected, the relative values 
y of the field at various points near a coil depend only on the shape of the coil. 
The value of the field at the centre of a coil is given by equation (4) and the 
» value at a point on the axis distant x from the centre by 


2rnl 
Ho, = 201(a,—a,) 
«| U+9) nat a, + {(i+%)? + a,?}1? (EIR a, + {(l— x)? + a 
Mn + (+a tay a, + (=a)? + a3) 


If (7) is divided by (4) and then both numerator and denominator of the right-hand 
side of the resulting expression are divided by a,, H),,/H 9 is obtained as a function 
of a,/a,, l/a, and x/a,. Thus, if in a design problem, instead of l/a, being 
specified, Hy,./Ho9 is specified for a given value of x/a, then, from equations (4) 
and (7), the value of //a, for each a,/a, can be calculated. However, the value 
of H,,,/Ho, is almost independent of a,/a, and thus it will often be possible to 
_ obtain the appropriate value of //a, for one value of a,/a,, and then use it with all 
required values of a,/a,. 

In general the magnetic field at points off the axis is not in the direction of the 
axis, but, with the rotationally symmetrical coils being considered, it will only 
have components parallel to the axis and in a radial direction. The general 
expressions for these components in terms of the value of the field at the axis are 


(Cosslett 1946) : 


ae he (= De (r\* 
: = —— (4) a (2n) 
Tn the direction of x : H,,=HAoz Z Ne hs +7 sip OFAN (5) ee 
See (8) 
F : Y if = 4/2 7 \ 2n-1 a 
and in the direction of r :H ,., = — 5 Hy, + Te ip Ss sta (5) Td eee) 
tee: (9) 


The differentiation is with regard to x. The coefficients can be obtained from 
equation (7) by differentiation but, in general, the full expressions for these and, 
consequently, the components of the field involve a very large number of terms. 
However, the case of the central plane of the coil («=(0) is much simpler. As 
the coils are symmetrical with regard to their central plane the value of the radial 
component of H,, is zero for all values of r. For small values of r the axial 
component may be represented by 


2 
r " 
— Ayo 


Ag =H — 4 
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where H, is given by equation (4), and, by repeated differentiation of equation (7) 
and putting x =0, 


oes 2nnl eat ait, Pe, RR, ae - 3 
00 ~ 10(ag— a) | (ap? + 2)? fa, + (aq? + 1?)¥?} (a,2+ 2)¥fa, + (a2 + 2?) 
i P 
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fe i 
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Note added in proof. Since the completion of this paper, a paper by G. Raoult, 
(1949) has come to the author’s attention. It is entitled ‘‘ Réalisation de champs 
magnétiques intense par impulsions. Applications aux phénoménes_ de 
polarisation rotatoire et de biréfringence magnétique”. In this, besides other 
work, Raoult describes coils, with the requisite theory, for producing half-wave 
intense magnetic field pulses. His theory for the design of coils is somewhat 
similar to that above but, in his final expressions, he has neglected the effect of 
the resistance of the coil on the peak current obtained. In other words, he has 
assumed that J=1. In Raoult’s work ky =(20a,//)¥7108S and D=1/10A"?. 
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ABSTRACT. A theory of contact noise is described in which the low-frequency noise is 
attributed to the random movement of adsorbed ions on the surface of a semiconductor from 
which an electron current is being drawn. Emission of electrons is assumed to take place 
only at localized patches on the surface and the adsorbed ions are assumed to give rise to a 
Schottky barrier layer, in which the potential maximum is linearly related to the concentration 
ofions. Diffusion of the ions over the surface gives rise to random fluctuations in the con- 
centration of ions in a patch which results in random fluctuations in the height of the potential 
barrier and the emission current. It is shown that for a circular patch the spectral power 
density of the noise current varies with mean current jg and frequency f as jo?f” over a 
small range of frequency and that x varies monotonically from —0-75 at the lowest frequencies 
to —1-125 at the highest frequencies. It is also shown that for a long thin rectangular patch 
the index x varies monotonically from —0-5 to —1-5 as the frequency is increased from zero. 
The dependence of the noise power density on temperature is also discussed. 


Se EN ERO DIV CALO IN 


N a previous paper (Macfarlane 1947) the author has described a theory 
of contact noise which depends on a diffusion conduction process in the 
atomic layers next to the emitting surface of a crystal. ‘The noise was 

attributed to the diffusion of clusters of mobile atoms on to the contact surface, 
their ionization on the surface and subsequent conduction away from the surface 
in the applied field. On the basis of a simple model it was shown that the spectral 
power density of the current noise should depend on mean current jy) and 
frequency f as j)”*1f-*. The index x, which lies between 1 and 2, depends on 
the concentration N of atoms in a mobile cluster. As x tends to 1, N tends to oo. 
However, experiments on silicon crystal rectifiers (Miller 1947), carbon granule 
microphones (Christensen and Pearson 1936, Miller 1949), germanium crystal 
rectifiers and lead sulphide photosensitive layers (Harris, Abson and Roberts 
1947) show that in some cases x is very nearly equal to unity and in other cases 
it must be taken less than unity. In addition some recently published measure- 
ments (Campbell! and Chipman 1949) on resistors have shown that the index 
of f varies from about —1 at low frequencies to — 1-6 at high frequencies. In 
order to explain these results on the above theory it would be necessary to 
assume enormous concentrations of atoms in mobile clusters. It is therefore 
apparent that this theory, in spite of its success in giving the observed dependence 
on current and frequency with micro-crystalline layers of PbS (Harris, Abson 
and Roberts 1947), is inadequate. 

In order to overcome these difficulties the author has developed a new 
theory, in which the noise is attributed to the random modulation of the 
Schottky barrier potential at emitting patches by random fluctuations in the 
concentration of mobile adsorbed ions. In its simplest form the theory refers 
to thermionic emission from a semiconductor, but it can be applied equally 
well to emission from one crystallite to another at contacts. 
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The idea is that the surface of a semiconducting crystal is covered by a partial 
film of adsorbed atoms, which have become negative ions by the capture of | 
electrons from the crystal. As a result a Schottky space charge layer is set up at 
the surface, as in Bardeen’s (1947) theory of surface states. The height of the — 
potential barrier is proportional to the concentration of ions. ‘The adsorbed ions 
are assumed to move at random over the surface of the crystal. When an electric 
field is applied to the crystal, emission of electrons is assumed to take place only | 
at localized patches on the surface. Fluctuations in the concentration of ions — 
in a patch occur due to the diffusion of ions out of and into the patch. These — 
fluctuations in number of ions in a patch give rise to fluctuations in the height of 
the potential barrier at the patch and therefore to fluctuations in emission of 
electrons. 

Since the electron emission is controlled by the concentration of ions in a 
patch the problem is to study the correlation in the number of ions in a patch and 
its dependence on time. 

The problem is essentially the same as a classical one in the study of the 
Brownian movement in colloidal solutions. It is to find the correlation function 
for the fluctuation in the number of particles in a cylindrical volume within a 
much larger volume of colloidal solution. An account of the relevant theory of 
colloidal statistics is given by Chandrasekhar (1943). 

In applying the formulae for contact noise to a multicrystalline layer or mass 
of semiconducting crystals it must be clearly realized that there are also other 
sources of noise, such as the Johnson noise of the bulk resistance, fluctuations 
in the number of current carriers, and noise due to the ‘ spreading’ resistance. 
The power density of these other sources is, however, independent of frequency 
at least at low frequencies, whereas contact noise per unit bandwidth increases 
as the frequency is decreased. If a frequency power spectrum of the total noise 
is known it is only necessary to subtract out the constant noise power, observed 
at high frequencies, in order to obtain the spectrum of the contact noise. 

In a very recent article, which came to the notice of the author after the 
theory to be described below had been worked out, Richardson (1950) has 
described in a most interesting manner a theory of contact noise in which the 
noise is also attributed to a diffusing adsorbed layer. He considers the region of 
multiple contact between two rough surfaces and assumes that the conductance 
per unit area is a function of the separation of the surfaces and the total 
concentration of adsorbate, near the point in question. By an analysis quite 
different from the one given below he derives an expression for the spectral 
power density which varies with frequency as f-!. Since the total noise power 
would diverge to infinity with this law Richardson suggests a modification to 
remove the divergence. With our model we arrive at a different expression 
which does not suffer from this defect. As we shall show it gives a law of the 
form jf” over a small frequency range, where the index w varies from —0-75 at 
very low frequencies to —1-125 at very high frequencies for a circular patch 
and from —0-5 to —1-5 for a long thin rectangular patch. 


§2. ANALYSIS OF THE PATCH MODEL 


We consider emission of electrons from a localized surface patch A of an 
emitter. Emission from the surface immediately surrounding A is assumed to 
be negligible. We assume that the current 7 through the patch A depends on the 
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ions can move at random over the entire surface of A and the surrounding surface. 
With these assumptions we can calculate the auto-correlation function and 
») from it, by Khintchine’s theorem, the spectral power density of the noise current 
‘jemitted by the patch in the following manner. 

The first step is to calculate the mean number of ions m in a patch A after 
‘time 7 from the instant when 7 ions occur in A. In doing this we assume that the 
/ motions of the ions are independent of each other and that all positions in A have 
)jequal a priori probability. Then if the probability that an ion somewhere in A 
will have emerged from it during the time 7 is P, the number of ions which on 
the average will leave the patch in 7 is mP and therefore the number remaining 
is n(1—P). In addition the number entering A on the average in + is NP 
» regardless of the initial value of n. The average number of ions in the patch after 
W time 7 is therefore 


Been oe Pe eee: (1) 


» Now in accordance with our assumptions the current 7 through the patch A 
depends on m as 7 =j,e“"-), where j, is the mean current. The noise current 


is therefore 
Maj, ie en im fpan— No 7 ee (2) 
) The approximation used in (2) is justified provided the noise current is small 


compared with the mean current. 
The auto-correlation function of the noise current is therefore 


f(r) =Aj(n).Aj(m) 
= (jga)*(n — N)(m—N) 


= (foa)*(n— N)*(1 — P) 

Vy ie Ot) te ner (3) 

In the third line of (3) we have averaged with respect to m keeping n constant 
and using (1), and in the fourth line we have used the well-known result that 

: the mean square deviation for random events is equal to the mean. 

: The next step is to consider the probability P. Suppose firstly that the patch 

- A is acircular disc of radius r and that the diffusion constant for the movement of 

/ ions on the surface is D, then it can be shown (Chandrasekhar 1943) that the 

* probability P that an ion somewhere within A will have emerged from it during 

the time 7 is 


(=e Whee ar (4) 
J 

where MCHA ai ae ake ae ene ER (5) 
and Le GO7 aay ORL | el eer (6) 


Using Khintchine’s theorem (Wang and Uhlenbeck 1945) we can now write 
_ down an expression for the spectral power density of the noise current. 


Res | ft) ARE SP ERAGNEAAC) A am anne (7) 
where F(p)= ip {1 — exp (— 4) [1,(a-¥*) + Ly(x-"4) |} cos(px)dx  ...... (8) 


and fe SOT, ES ENS Ee ee ee rr eee ee) 
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is shown in the Figure. 'T'wo limiting formulae are useful. 
When p<K<!l E(p)=20'2345p7 25 = = eee (10) 
and when p>! FO\=O 1406p" ec eee (11) 


In the preceding analysis R(w) is the mean square noise current at a single § 
localized patch or contact. It was derived on the assumption of a constant voltage | 
across the contact. When there are many such contacts, as in a multi-crystalline : 
resistor or at contacts between rough surfaces, the noise fluctuation of the total | 
input current is a weighted sum of the separate fluctuations at the contacts. 
The weighting factor will depend on the current at the contact and on the current 


10 


Long Thin Patch 


Circular 


103 107? 107! | 10 10? 
P= WT. 


which a series E.M.F. at the contact would produce at the input to the resistor. 
It will therefore depend on the position of the contact in the resistor. 

As an example consider a resistor which can be represented by m sets of © 
contacts in parallel, each set consisting of m contacts in series. Between adjacent - 
series contacts assume a ‘ bulk’ resistance. Let the resistance of the pgth contact | 
be r,,, and the ‘ bulk ’ resistance in series be R,,,. ‘Then a simple circuit analysis 
shows that the mean square noise current at input is 


n 


Ape =jeR? x S Coif fiR,. Me ae (12) 
p=l1q 


pad 
=1 


where jy is the mean current at input, R the total resistance of the resistor, 
n ‘ 


R,= & (Ryqgt%pq)and 
q=1 


: : 


bog HAG NG Tog! (Oing) a. » ae (13) 


pa pa 


As a further simplification suppose that all the contacts have the same 
resistance 7, and all the resistances R,,, are equal to Ry. Then 


Aj? =jo*{ro/(Ro +70)}*e/mn. aha 
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Although this last expression applies only to a very simple resistor it does indicate 
that the power density of the noise current decreases as the number of contacts 
‘is increased. 

) On account of the summation in (12) the variation of noise power density 
jwith frequency will not be quite given by (8). However, at high enough 
mrequencies where (11) is valid, c,,~0-5864a?N tp /®w-*® and Aj? is then 
‘directly proportional to w~®8, Similarly at low enough frequencies Aj? varies 
‘with frequency as w-*/, 


Seas lH Pp OEPRECT OF PATCH SHAPE ON THE SPECTRUM 


In §2 we assumed a circular patch. However the shape of a patch does have 
an interesting effect on the spectrum, as we can see by evaluating R(w) for a long 
‘thin strip. Then diffusion takes place primarily across the strip. If the width 
' of the strip is / the probability P is given by (Chandrasekhar 1943) 


1 — P=erf («-¥2)—(sx/m)Y[1—e-V*], saa (15) 


Sewhere x=7,/7 and 7,=h?/4D. 
The spectral power density of the noise current is then 


RO) = 4 aaNet eee PMU esters (16) 
here p=wrt, and 


i p)= ie ferf (x—1?) — (x/a)"?[1 — e-V7]} cos(px)dx. ~...... (17) 


The function F(p) has been evaluated from the formulae 


Zz co (—)"a/rp” 3? 

| NP)= = Wat il a+iysinaGn— 
| zs < (—)?tha/mp?” (18) 
/ 


y—o 21 (2n + 2)P'(2n+ 5/2) 
Phy Webtort pete (19) 
. 
sand open lor peels oye 2 ES (20) 


F(p) for a long thin strip is shown plotted in the Figure for comparison with F(p) 


| These formulae were derived by the method given for (8) in the Appendix. 
} 
| for a circular patch. 


§4. DISCUSSION OF THE NOISE SPECTRUM 


Consider firstly the spectrum for a circular patch. If F(p) is replaced by a 
power law cp* over a small range of p then it can be seen from formulae (10) 
/and (11) and from the Figure that « varies from —0-75 when p<1 to —1:125 

when p>1. From this it follows that the integral of F(p) taken over the entire 
frequency range from p=0 to p= © is finite, notwithstanding the fact that F(p) 
tends to infinity as p tends to zero. 
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It will be observed that the spectral power density of the noise current 
increases with the square of the mean current and directly with the average 
number of adsorbed ions in a patch. 

Only an approximate expression for the dependence of R(w) on temperature} 
can be given unless the variation with temperature of the mean number, N, of 
ions in a patch is known. If we assume N independent of temperature T we have 


a=EjkT, D=D eH, — r9=(2r)?DeoteB*?, 2.2... (21) 


where E is the activation energy for diffusion of ions over the surface and KE, 
is the change in height of the barrier potential at a patch produced by one ion. 
Then, if wr) >1, we get from (7) and (11) that 
R(w) ~0-186877,2ND,!8(2r)-14(Eg/RT Pe oh 98 
oo te eRe ee (22) 


From (22) it follows that R(w) increases with T for temperatures below| 
T,=E/16k and decreases with increase of T for temperatures above T,.. 
If E=0-lev., which is probably rather a high value for surface diffusion, | 
T,=60°x. Then for temperatures above 60°k. and for high frequencies R(w) 
should decrease as the temperature is increased. Evidence that R(w) does 
decrease with increase of temperature is provided by measurements on silicon 
rectifiers (Miller 1949). 

Although little is known of the values of the diffusion constants D, and £, 
for diffusion of adsorbed ions over a crystal face both Dy and E are likely to be ! 
much smaller than for diffusion of impurities through a crystal (Mott and Gurney * 
1948). In order to appreciate the magnitude of 7, take D)=10-*cm?/sec., , 
E=0-lev., and take the radius of a patch as 10-*cm. Then 7=0-016sec. at | 
room temperature. 

Consider now the effect of patch shape on the spectrum. We have found that, , 
when wt)>1, R(w)xcw-*® for a circular patch and R(w) «<w-*? for a long thin | 
rectangular patch and that, when wr)<1, R(w)«w** for a circular patch and | 
R(w) cw-\? for a long thin rectangular patch. Deviations in the index x in the : 
expression R(w)ocw* can therefore be explained in terms of patch shape so long. 
as x lies in the range —1-125>x>-—1-5 at high frequencies and in the range 
—0:75<x<—0-5 at the lowest frequencies. It has been found experimentally — 
(Miller 1949) that for a photoelectric cell illuminated with a constant intensity 
of light «= —1, for a carbon granule microphone — 1>x> — 1-1, and for a carbon 
resistor —1:1>x>-—1-2. These results can therefore be explained by assuming 
nearly circular patches. For various types of resistor Campbell and Chipman 
(1949) have found that the index x varies from about —1 at a frequency of 
20 kc/s. to about —1-6 at 500kc/s. and that R(w) increases as the square of the 
mean current. It would therefore be necessary to assume long thin patches in 
order to explain these results. For micro-crystalline layers of PbS it is found 
(Harris, Abson and Roberts 1947) that x may lie between —1 and — 1-4. 

It appears therefore that the theory of contact noise described above, in which 
the noise is attributed to fluctuations in emission from localized patches on the 
surface of an emitter due to fluctuations in the number of ions in a patch, can 
provide a reasonable explanation of the observed results. 
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APPENDIX 
EVALUATION OF THE INTEGRAL FOR K(p) 


The Mellin transform of 1 —P(x)=1— exp(—«-™4)[Iq(a- 4) + 1,(«-¥4)] is 
(Macfarlane 1949) 


= 28-21 ( + 45)0(—45) 


0% \/nl'(2 +45) 


0<o<h. 
By using the inversion formula of the Mellin transform to express 1—P as an 
» integral in the s-plane and reversing the order of integration we get 


: [1— P(x)] cos (px)de = 5 fe  taeee fies Tee 


_] etic 1 oy 
ALTON My coc 
i The integrand has poles at 
s=—4(n+4)= —-}, —% —3,.... due to [($+4s) 
and at s=n/4 (n>0) =},4,3,..-. due to I'(—4s) 
| and double poles a 
Ret eee 1, 3, Dees oe due to I'(—4s) and ['(4 — $s). 


: Moving the contour steadily to the right we can take out the residues at the 
poles to the right of o and get the following convergent series for F(p). 
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ee Det +2) G2) cos (a8) 
n~8r+4 
5 Ae ea eet 
n—o V7 (8+ 6) (8" + 5)P(2n+ 1) —Y(2n + 1) +1n(16p) 


The first few terms are 


0-23447p-84 — 0:31333p-" + 0-34892p-¥4 — 0-14378 +0-03646 In p 
— 0-04953p"4 + 0-00718p"2 — 0:00123p34 + . 


: Moving the contour to the left we obtain an asymptotic series for /(p) from the 
) residues at the poles of ($+ 4s). 


4)P(n—4)T(4n +8) 


(—)0@+ 8) PRs et 
Kpy~ >, 721 (n+ 1) b sin. (ee) 2 ee pene 


~0-14659p-98 — 0-04925p-18 — 0-01394p-138 — 0-01093p-158 
—0-01328p-178 ... 
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Dielectric Loss and Dielectric Constant Measurements in. 
Supercooled Liquids é 


By C. DODD ‘anpb: G.-N. [ROBER ES 
University College and Queen Mary College, London 


Communicated by E. N. da C. Andrade; MS. received 26th Fune 1950 


ABSTRACT. In order to provide further evidence in support of the structural change © 
which takes place in a liquid when it is supercooled, measurements of dielectric loss at high | 
frequencies and of dielectric constant have been made on various polar liquids. The ~ 
evidence from the dielectric loss measurements is inconclusive but for the four liquids tested _ 
a significant discontinuity in the dielectric constant—temperature curve has been established — 
at the melting point in each case. 


SiN RO DUCEION 
HE work of Dodd and Hu Pak Mi (1949) has shown for certain polar liquids 
| that, although there is no discontinuity of the viscosity 7 at the melting point, 
there is a sharp change in the course of the line of log y plotted against 1/T | 
at this point. Such a phenomenon can only result from some sort of structural 
change in the liquid as it supercools and any such change should also show itself 
in the temperature variation of other properties of the liquid. 
The present paper describes measurements made on the dielectric loss and 
dielectric constant of several liquids over a range of temperature on both sides of 
the melting point. 


$2. DIRE Cap ICS WOS)s 
‘The apparatus for the determination of the dielectric loss was one in use in the 
research department of Electrical Engineering at Queen Mary College, London. 
Oscillations of wavelength 3 cm. from a klystron were fed into an Hy, cavity | 
resonator containing the liquid, and were passed through coaxial cable to a crystal 
detector, variable attenuator and galvanometer. The resonator was in a constant 
temperature bath and a thermocouple in contact with the movable piston attached 
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» to a micrometer shaft registered the temperature of the liquid. The width of the 
“sresonance curve was obtained from the positions of the micrometer shaft when the 
y galvanometer deflections were half the maximum value. From this width the 
dielectric loss was calculated. In Figure 1 the variation of loss with temperature 
~ is shown by plotting e«”/e) against temperature for both normal and supercooled 
liquid phenyl ether (M.Pt. 26-9°c.) where e” is the imaginary part of the complex 
‘) dielectric constant and e, the dielectric constant of free space. It is seen that over 
“this range the scatter of the points is too great for any indication of a discontinuity 
+ at the melting point to be apparent. 
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Soe DIRLEC TRIG CONSTANT 

In order to measure small changes in the dielectric constant of a liquid to a 
| high degree of accuracy, a test condenser was joined in parallel with a standard 
) and a variable condenser in a circuit tuned to a frequency of 244 ke/s. 

The test condenser (Figure 2) was of the parallel plate type, the plates P,, P, 
) of diameter about 6 cm. being separated by a distance of the order of 1 mm. by the 
; quartz washers QQ. ‘The empty capacity was determined at various temperatures 
» by amodified Hartshorn circuit. ‘This value was about 28 pr. and changed by only 
0-18 pr. for a temperature change of 40°c. The capacity varied regularly with 
» temperature and showed no hysteresis when taken over a heating cycle. The 
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whole cell was then filled with liquid, the thermometer T registering the temper- . 
ature to an accuracy within 0-02°c., and was immersed in a bath whose ; 


temperature could be maintained constant to better than 0-01°c. 


The variable condenser was of the cylindrical type, the capacity being varied _ 
by movement of the inner cylinder. The total movement available was 1,500 
divisions on the micrometer, one division (10-3 cm.) corresponding to a capacity © 
change of 0-01074 pr. A change of half a division caused a measurable change in | 


the resonant position. , 

The setting of the micrometer condenser to give resonance was found with the 
test condenser at various temperatures both above and below the melting point 
of the liquid. The liquids chosen, which were all polar, readily supercooled and 


remained in the liquid state throughout. The micrometer condenser setting C is _ 
plotted against temperature for phenyl ether, M.Pt. 26:9°c., w=1-17 (Figure 3); | 
salol, M.Pt.42°c., 4 =3-15 (Figure 4);-menthol, M.Pt.42°c., »=1-56 (Figure 5)m) 
and azoxybenzene, M.Pt. 36°c., »=1-70 (Figure 6); yw is the dipole moment, in — 
all cases x 1018. Because of the large ionic loss for azoxybenzene at the lower 


frequency (244kc/s.), the circuit was tuned to a frequency of 1-65 Mc/s. for this 


liquid. Values for the dielectric constant K deduced from the corresponding — 


micrometer condenser settings, are given as ordinates on the right-hand side of the 
graphs. 
§4. DISCUSSION 


Examination of Figures 3, 4, 5 and 6 shows in each case that for the liquid above 
its melting point the points lie fairly well on a straight line. Those for the super- 


Melting PHENYL ETHER 
Point. 


300 |- 


Deviation 
8 


Melting 
% Point 


200 |- %s 


10 20 30 
Temperature (°C,) 


Figure 3. 


cooled liquid also lie on a line which is, however, inclined to the line for the norma 
liquid. In each case these lines intersect in the neighbourhood of the melting 
point. ‘This discontinuity at the melting point is brought out more strikingly by 
the insets to the graphs, where the deviations of all points from the straight line, 
fitted by least squares, for the liquid above the melting point, are plotted against 
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temperature. It should be mentioned that in order to avoid any systematic error 
arising in the apparatus, the dielectric constant was measured first in the region 


above the melting point, then in the slightly supercooled region, back into the 
normal region, further into the supercooled region, and so on. 
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For phenyl ether, salol and menthol similar discontinuities at their melting 
_ points have been found by Dodd and Hu (1949) for the variation of viscosity with 


.temperature. ‘The experimental technique used for these viscosity measurements 
| 
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is not, however, easily adaptable to azoxybenzene, which is a dark, almost opaque 
liquid. 
It would therefore appear that the change occurring in the nature of the liquid’ 
as it enters the supercooled region is responsible for both the discontinuity in 
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dielectric constant and the discontinuity in viscosity. ‘The search for similar 

discontinuities in other physical properties of these supercooled liquids 1s 
continuing. 
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The Dielectric Strength of Aluminium Oxide Films 


Frohlich (1937) has shown that, under the influence of an applied electric field F,’ 
electrons in a crystal will gain energy A from the field and lose energy B to the lattice ‘ 
vibrations, and that the resultant energy change B—A is given by B—A=B—(e?/m)rF”,’ 
where e and m are the electronic charge and mass respectively, and 7 is the average time} 
between collisions. Electrical breakdown occurs when electrons possessing energy greater} 
than the ionization energy J gain energy from the field. Accordingly, the breakdown field 
strength is inversely proportional to the square root of the relaxation time. Hence Fa/-??,/ 
where / is the mean free path of the electrons. If the latter is decreased by introducing} 
foreign ions into the crystal, or by raising the temperature, the dielectric strength increases. | 
In a thin film the value of B increases as the film thickness approaches the mean free path, 
since the number of collisions with the film boundary increases, and thus the dielectric} 
strength of thin films of dielectric should rise with decreasing film thickness. Fréhlich has! 
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calculated the order of such an increase for mica. Austen and Whitehead (1940) using mica 
) of thickness varying from 10~* to 10~° cm. have verified his predictions, and more recently 
Plessner (1948) working with CaF, NaF and KBr films has also confirmed Fréhlich’s theory. 

This theory may be extended to films of amorphous materials, such as aluminium oxide, 
for which there should be a similar increase of strength although, since the electronic mean 
free path in amorphous substances is considerably less than in crystalline material, the effect 
should only occur to any marked extent in very thin films, for example films of thickness. 
probably less than 100 a. 

The aluminium oxide films used in the present investigation were made by partially 
anodizing a thick aluminium film, formed by evaporation of the metal on an optical flat in a 
y yacuum of 10-*mm. Hg. The metal film was partially immersed in an aqueous solution of 

ammonium phosphate, and to prevent ‘ peeling’ of the film the flat was subjected to ionic 
+) bombardment before evaporation. 
The thickness of the oxide film was measured by employing equal chromatic order 
interference fringes. A cross section of the film is shown in Figure 1. It projects beyond 
+ the adjacent aluminium film, so that the oxide film may be defined by the two distances 
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} Figure 1. Figure 2. 


| tits thickness and t, the step between surfaces A and B. Hass (1949) has shown recently 
) that the ratio t,/t is constant, independent of t and equal to 0-27. Thus the optical path 
’ difference between rays (i) and (11) reflected from the air—aluminium and the oxide— 
» aluminium surfaces, respectively, is given by 2(u—0-27)t. Taking w=1-:65—as given by 
+ Hass—the optical path difference between the rays is 2:76t. When placed in contact with 
a semi-silvered optical flat, the step may be detected by means of fringes of equal chromatic 
» order, formed by reflection at the film-optical flat surfaces, white light and a constant 
» deviation spectrometer being used to produce these fringes. Since a thin surface oxide 
layer exists at the air—aluminium interface, phase changes at the oxide-aluminium and 
air—aluminium interfaces are equal. ‘The fringes on both sides of the step were parallel and, 
/ accordingly, the oxide films were of uniform thickness. 
Attempts were made to measure the film thickness directly by completely anodizing 
) aluminium deposited on a film base of gold, nickel or platinum. With this arrangement the 
i film thickness could be measured independently of the the refractive index and the constant 


¢ 


~a 


| t,/t but such metals ‘ peeled off’ the flat during electrolysis. 
The potential difference across the oxide film was applied by means of two spherical 
| electrodes—steel ball bearings + z,q in. in diameter—one in contact with the surface A and 
the other lightly sprung so as to rest on B. ‘The pre-breakdown currents, measured with 
: gradually i increasing potentials, were similar to those obtained by von Hippel (1938) with 
, mica and alkali halide crystals. "They are probably produced by strong field emission from 
, the cathode although it has been suggested that they may be formed by ejection of electrons 
. from impurity centres at high field strengths. 
At a certain voltage this pre-breakdown current suddenly increased to a very high value 
and this voltage represents the breakdown value. Many measurements were made with 
the same film at different places on the film and averaged. In all experiments the film was 
in a vacuum and at 15° c. Figure 2 shows the variation of dielectric strength with film 
thickness. As the latter increased from 100 A. to 1,000 a. the strength decreased from 
11 x 108 v/em. to 610% v/em. The value at 1,000 a. agrees with that found by Guiather- 
schulze and Betz (1936) for aluminium oxide in bulk. 
3 K-2 
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The dielectric strength of one film was measured at 15° c. and at 100° K. It increased 
from 8108 v/cem. to 9°510%v/cem. Similar results were found by von Hippel and 
Maurer (1941) using soda lime glass. 


amorphous structure and supports Hass’ statement that no pores exist in the oxide film 
formed by anodization. 


The relationship between formation voltage and film thickness in the present experi- 


ments is similar to that found by Hass (1949) and by Giintherschulze and Betz (1936). It _ 
is linear, but whereas the latter obtained a slope of 11 a/volt and Hass 13 a/volt, the preseges 
work gives 14 a/volt. 
The Physical Laboratories, P. D. Lomer. 
The Washington Singer Laboratories, 
University College, Exeter. 
7th July 1950. 
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Ferromagnetic Resonance in Manganese Arsenide 


Manganese arsenide is known to have a ferromagnetic critical point in the region of 
40° c. and to exhibit pronounced temperature hysteresis (Bates 1928, 1929, 1931, 1933). 


Below the critical point ferromagnetic properties are clearly marked, and at room temperature : 


the material has a saturation intensity of magnetization comparable with that in nickel at 
the same temperature. 
absorption in this material at a wavelength of 1:25 cm. and at temperatures from 18 to 50° ¢ 


The intensity of magnetization in small rod samples was measured by a simple ballistic : 


method throughout the required temperature range and in applied magnetic fields up to 
5,000 oersteds. It will be seen that at room temperature a value of about 400 £.M.U. is 
found for J, whereas Bates (1933) obtained a value of about 600; in each case satisfactory 
checks against a pure nickel sample were made. The reason for this discrepancy is being 
investigated. 


This indicates that aluminium oxide has an 


An investigation has been made into the ferromagnetic resonance * 


emis 


; 


From each rod, very thin discs were cut and mounted to form part of the lower end ofa . 


cylindrical cavity resonating in an H,,;, mode. The power absorbed in the specimen, at 
constant microwave frequency, was measured as a function of the strength of an applied 


steady magnetic field. The applied field, Hyax, for maximum power absorption was thus | 


determined at each temperature. Using the relations first developed by Kittel (1947, 
1948), it follows for a specimen of the form used that 

2X 10°7 
Aa( BA )max*!* : 


Q 


Z= 


where g is the Landé splitting factor, A, is the free space wavelength in centimetres, and 
B and H are expressed in kilo-oersteds. 

Specimen results obtained with one sample are shown in the Table. Near the critical 
point the intensity of magnetization is changing very rapidly with temperature, the total 
absorption is small, and the experimental accuracy is therefore decreased. 


Ne 255 cm. 
t (C5) 200/298 533"1) 235-08 936-00 S7-0Mes7 28ers 9 Olms 9a7S 41:7 45 and 50 
FAynax (0e-) 2700 2725 2800 2800 2925 3050 3280 3650 4150 ~4550 4500 to 5000 
{ (8.M.U.) CW seeeoel) koe Sisk. So SLOP IAG) 50 = — 
g SS 3°75 Seda S75 5:65 SnO¢ 3°65 5-0F 38, S07/G Se to 3°85 
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Three different specimens were examined and, within the limits of experimental error, 
e g-value determined by the above equation does not vary up to the critical point. At 
room temperature, the three g-values found were 3:7., 3°7;, and 3°85. 

It is of interest that resonance absorption could be detected above 45° c. where the man- 
‘}ganese arsenide has paramagnetic properties. 'The absorption is small and broad, so that 
the position of the peak can only be found approximately, but it appears to be such as to 
suggest little, if any, change from the ferromagnetic g-value. There is no discontinuity 
»tobservable at the critical point, the width of the absorption curve slowly increasing with the 
temperature. At 50° c., the width is about 20°% greater than at room temperature. These 
» latter observations are in accordance with those of Bloembergen (1950) whose results for 
») nickel and supermalloy were published during the course of this investigation. 

> Of particular interest, however, is the very large g-value which is obtained. Previously 
~) published results by various authors (see Kittel 1949) have given values mainly of the 
)) order of 2:20 and the excess of this above the value 2-00 for an electron spin has been 
5 qualitatively interpreted in terms of the contribution of the orbital moment (Kittel 1949, 
Polder 1949, Van Vleck 1950). Kittel has shown that if g is the g-value deduced from 
~) microwave absorption measurements and g’ that from direct gyromagnetic measurements, 

» g should be greater than 2-00 by the amount that g’ is less than 2:00 for small variations 
about the electron spin value. 

If this type of explanation is to hold in the case of manganese arsenide, one would expect 
g to be very much less than 2:00. It is thought that there has been no such measurement 
on manganese arsenide itself, but Galavics (1939) has reported measurements on an Mn—Sb 
} alloy stated to be Mn,Sb, for which he obtained a value for g’ between 1-9 and 2-0. 
It is hoped to continue this investigation with other materials of the types Mn—P, Mn-Sb, 

1? Mn-Bi, and results will be reported later. 

The authors wish to record their thanks to Professor L. F. Bates, who kindly supplied 
+ the manganese arsenide samples, for his continued interest and advice. 
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25th July 1950. 
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Electrical Conductivity of Gray Tin 


| Gray tin has a non-metallic crystal structure similar to diamond (Bil 1919, Brownlee 
+ 1950), i.e. it is the end-member of the series carbon-silicon-germanium—gray tin. One 
} would therefore expect its electrical properties to be appropriate to its position in the series ; 
/ but there seem to be no reliable experimental results available in the literature. It was 
/ originally reported that gray tin was a superconductor, but de Haas (1928) and Sharvin 
| (1945) showed that this was incorrect, at least above 1:2°k. Busch (private communication) 
has obtained some values for electrical conductivity of gray tin powder which show it to be 
- asemiconductor. On the other hand Moesveld (1937) reported that it was a normal metallic 
' conductor. It seemed unlikely that results of real value could be obtained until a way was 
found of preparing gray tin in massive crystalline form free from white tin. Normally 
it is prepared by keeping tin at a temperature well below its transition point (13:2° c.) 
when the metal crumbles into a friable powder. The large change in specific gravity 
from 7:31 to 5-75 is no doubt the cause of this disintegration. Unfortunately, subsequent 
compression of the powder tends to ‘ransform it back to white tin, even when the com- 
P-ession is carried out at low temperature. 


——— 
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In transforming a kilogram of spectroscopically pure tin it was found that amongst - 


the large quantity of fine powder produced were a few lumps of gray tin several millimetres 
in linear dimensions. Suitable electrodes were pressed against these lumps and their 
electrical conductivity measured for small D.c. currents, using a milliammeter and potentio- 
meter. No contact non-linearity was observed, but precautions were necessary to eliminate 
the rather large thermo-electric E.M.F. generated under non-uniform temperature conditions, — 
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a-tin: variation of conductivity with temperature. 


The variation of conductivity with temperature is shown in the Figure; the values given for 
specific conductivity must be regarded as only approximate as there was some difficulty in 
measuring the cross section and electrode areas of the specimens. Almost identical results 
were obtained for three different specimens. 

It is seen that gray tin is an intrinsic semiconductor, in which the conductivity is given by 


o=A exp (—€/2kT),. 


From the upper portion of the curve we obtain d~6 x 10* ohm cm“, and e=0-098 ev. 
At the lower temperatures the presence of impurities in the tin is presumably causing 
enhanced conductivity and a lower activation energy. Gray tin thus fits consistently into the 
series of intrinsic homopolar semiconductors as shown in the Table below. The very small 
energy difference between the filled and unfilled bands is noteworthy, and leads to the high 
conductivity at 0° c., remarkable for a semiconductor, of approximately 2°5 x 10° ohm7! cm. 


€ A 

Silicon 1:12 0-9 x 102 

Germanium 0:75 Boo Oe 

Gray tin 0-1 6 x 104 
The author wishes to acknowledge the award of a Leverhulme Research Fellowship. 
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Development of the Low Pressure Electrodeless Discharge 
in a High-Frequency Electric Field 


The starting of a high-frequency electrodeless discharge with the mean free paths 
of electrons much larger than the dimensions of the glass vessel has recently been 
3 investigated (Gill and von Engel 1947, 1948). In these experiments a cylindrical glass 
»/ tube containing gas at this low pressure was placed in a uniform high-frequency electric 
> field (wavelength order 10 m.) directed parallel to the axis. It was shown that, if an 
electron starting at one end wall with a small energy and in a suitable phase of the field 
) crossed the tube in one half-cycle, and hit the opposite end wall fast enough to release 
7) more than one secondary electron, then an electron avalanche would develop. The 
+ starting field strength should therefore depend on the material of the end walls, and not 
fon the gas. This was confirmed by experiment. As the discharge develops a current 
passes through the gas and the field decreases. This working field strength, and the 
spectrum emitted, are characteristic of the gas, and do not depend on the material of the 
~ vessel. 
It was derived that in the initial stages of the discharge electrons have maximum 
“> energy (about 100 ev.) when they hit an end wall. This is enough to release between 
» one and two secondary electrons per incident electron, most of which contribute to the 
+» development of the avalanche, whilst a few are lost to the side walls. This maximum 
energy is also much greater than the ionization energy of any gas. 
But in the fully developed discharge both photographic measurements at wavelengths 
of about 3,000a. and visual observation show that most of the light is emitted from the centre 
of the cylindrical vessel and very little from the ends. Assuming that this light comes 
directly from atoms excited by collisions with electrons, it appears that the electrons are 
fastest in the central region. A cloud of electrons seems here to oscillate with small 
amplitude about a space charge of positive ions in the centre of the tube. Very few of 
these electrons hit the end wall; if they did they would be lost. But it is known that 
when losses to the walls are large (e.g. beyond cut-off wavelength (Gill and von Engel 
1948) no field, however large, can start the discharge. 
As far as we know no atternpt has been made to describe the transition from the initial 
to the final stage of such a discharge. At first sight it would appear that, of the few ions 
formed in the gas in the earlier stages, most would be produced near the walls where the 
electrons are fastest. ‘Two clouds of positive space charges would then be formed near 
the end walls, and this would continue to make the electrons fastest near the ends of the 
tube. (The effect of the positive wall charge caused by secondary electron emission 
» has been considered, but found irrelevant here.) 
f This, however, is not true. A calculation of the energy of an electron as a function 
of position along the tube shows that any electron capable of starting the discharge 
acquires energy equal to the ionization energy of any gas well before it reaches the centre 
» of the tube. The ionization produced is obtained by combining this result with known 
( ionization probabilities. 

Now electrons crossing the tube in one direction produce a distribution of ions in 
space which increases rapidly towards the centre, and then slowly from the centre to the 
) far wall. Electrons crossing in the opposite direction produce a symmetrically opposite 
| distribution. The sum gives the maximum number of ions at the centre. As the 
¥ rnultiplication process continues the number of ions increases, but the maximum 
) ionization is always at the centre. 

When the number of ions becomes considerable the field in which electrons move is 
| the applied field plus the field due to this positive space charge. As the latter increases 
_ the electrons reach their maximum speed progressively nearer the centre of the tube. 
_ New electrons are now produced (by ionization of gas atoms) near the centre, and start 
oscillating in a phase different from that of the main cloud. Some are lost to the walls; 
most oscillate in the gas with a maximum speed near the centre, their amplitudes becoming 
- smaller the larger the positive ion space charge. These are the electrons which are found 
in the final stage of the discharge. The electrons originally produced by secondary 
emission from the walls are progressively lost to the end walls as their speed at the ends 
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of the tube is reduced by the influence of the positive ion space charge. ‘These electrons ; 
neutralize the positive wall charge. i 

Positive ions in the gas drift to the walls owing to self-repulsion in the ion cloud. 
However, electrons produced in the gas in unfavourable phases of the applied field are ; 
also lost to the wall, neutralizing these positive ions and charging the walls ligne 
negatively. 

Thus it appears that the transition from multiplication by secondary emission to an © 
equilibrium state sustained only by ionization in the gas can only occur when the electrons | 
produced by the first process are all replaced by electrons produced in the gas which — 
then oscillate with an amplitude smaller than the length of the vessel. 


The Clarendon Laboratory, G. FRANCIS. 


Oxford. A. VON ENGEL. 
10th August 1950. 


GILL, E. W. B., and von ENGEL, A., 1947, Nature, Lond., 159, 404; 1948, Proc. Roy. Soc. A, 192, 446. 


REVIEWS OF BOOKS 


An Introduction to Heat Transfer, by M. FIsHENDEN and O. A. SAUNDERS. 
Pp. x+205. 1st Edition. (Oxford: The University Press, 1950.) 15s. 


There seems to be a general trend towards putting the engineering treatment of heat~ 
transfer problems on to a more systematic basis than has been the practice in the past, and 
the authors of this book have been among those who were instrumental in bringing it about. 

They have now published a book, distinct from their earlier joint writings, in which the 
whole subject is surveyed on a reasoned basis. 

The first chapter deals with radiation, whether the theoretical simple radiation from a 
black body or the complicated phenomena encountered, for example, when considering the — 
radiation from the flames in a furnace. The questions which arise when geometrical 
considerations intervene, or when the emitting and receiving surfaces are of different 
emissivities, are considered. 

The next chapter is on conduction, and does not deal in detail with the calculation of 
heat flow in complex systems, though it does treat the theory of sufficient cases for most 
practical applications. ‘The theory is a matter for mathematicians, and is treated in separate 
books, of which that by Carslaw and Jaeger is probably the best known, though a recent one 
by Ingersoll, Zobel and Ingersoll should not be overlooked, especially as it and the book 
under review are in a sense complementary, and share the same outlook—that of keeping 
basic physics in view while applying the theory in practice. In this chapter, Fishenden and 
Saunders describe in detail the graphical method due to Schmid for calculating heat flow, 
and they also refer to the use of relaxation methods, and to electrical analogies. 

About 60% of the book is devoted to convection, first the general theory, then the relation 
between convection and friction, then forced and natural convection in the commonest 
shapes of system, concluding with convection to and from fluids which evaporate or condense 
in the course of the action. 

The material is admirably presented, and is illustrated by worked examples, which add 
considerably to the value of the exposition. Many useful charts and tables are given, which 
will eas» the work of calculation in specific applications to real systems. J. HA. 


Die Ubermikroskopie, by B. von Borrirs. Pp. 416, with 70 plates and 225 line 
diagrams. 1st Edition. (Aulendorf, Wiirttemberg: Editio Cantor, 1949). 
48 D. Marks. 


For some ten years now there have been available commercial models of the electron 
microscope which can far surpass the optical microscope in resolving power and which, 
although several times as expensive as the latter, are not exorbitantly priced as modern 
physical apparatus goes. Nevertheless the amount of research work of fundamental 
importance with which it can be credited remains relatively small. After considerable 
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~) enthusiasm, in principle, over the introduction of the new instrument, there has followed 


in many quarters a certain caution in regard to its practical application. This degree of 
disillusionment has been due in some part to a lack of understanding of the natural limitations 
attaching to its technique, so that an immediate extrapolation was often expected to be 
possible from optical practice into the smaller orders of dimensions. In fact the preparative 
methods, and the way in which image contrast arises, are so different in electron microscopy 
that quite new interpretative methods are required. This is not to say that a great deal of 
important progress has not already been made in metallurgy, biology and most of the applied 
sciences with the aid of the electron microscope. But much more might have been possible 
if a clearer realization had existed of the natural limitations to its use. 

It is with the exploration of these limitations that Dr. von Borries, who was responsible 
with Dr. E. Ruska for the original development of the magnetic electron microscope, has 
principally concerned himself in recent years. An account of these researches forms the 
core of his present book. He has investigated the instrumental factors of illumination 
brightness, which determines the maximum magnification which can be used and hence the 
resolution limit, the recording sensitivity, which sets a lower limit to the total amount of 
energy passing through the specimen during the exposure necessary to form an image, and 
the fluorescent efficiency of viewing screens, which determines the lowest permissible 
illumination for focusing and hence the temperature rise in the specimen. In the light of 
these findings he examines the possibility of electron microscopy of living material, as 
limited by the thermal and ionizing effects of the electron beam. The lethal dosage is so 
quickly reached that observation would only be practicable at low magnification, and hence 
at low resolving power; only in the case of a few organisms could the optical limit be 
exceeded, even if a sufficiently high voltage were available to penetrate them. 

Aside from the biological aspect of the enquiry, with its pessimistic conclusion, the 
author presents a great deal of valuable information as to the efficiency of electron guns, the 
scattering and energy loss which accompany the passage of electrons through matter, and 
the characteristics of photographic emulsions towards electron illumination. ‘The present 
limit set to resolving power by the aberrations of electrons is also fully discussed, although 
possible means of correction are hardly mentioned. 

As the book is primarily intended for scientists who may wish to assess the positive value 
of the electron microscope in their particular field of research, the detailed exposition of its 
limitations is preceded by a review of its principles and constructional features, and is 
followed by a survey of its applications in pure and applied science, the three sections being 
of approximately equal length. An appendix has been added to bring the treatment up to 
date to the beginning of 1949, the main text having been completed a year earlier. ‘The 
result admirably fulfils the author’s intention of giving a well-balanced account of the 
present state of development of electron microscopy. Omissions of reference to non- 
German work, which in the difficult war and post-war conditions of communication might 
have been expected to be serious, are relatively minor, although it is admitted that in many 
cases the author has not had access to the full text of a paper cited. he illustrations are 
naturally almost entirely from German sources, and thus permit a much better assessment 
of the progress made in electron microscopy there since 1939 than was possible from the 
compressed and unillustrated FIAT reports. 

The quality of the reproduction is high, but the same cannot be said of the rest of the 
technical production of the book. It is probably inevitable that the quality of the paper for 
the text should be poor, in present-day circumstances, but the binding could surely have 
been better executed. The pages have been inserted individually, not in stitched and folded 
sections, so that they come loose with very slight use. It is regrettable to have to make 


such a complaint as to the production of what is in all other respects an excellent monograph. 
V. E. COSSLETT. 


Colours and How We See Them, by H. Hartrrince. (Royal Institution Christmas 
Lectures 1946-47). Pp. xi+158. 1st Edition. (London: G. Bell & Sons 
ed 1949). 15s: 


This book is based on the Royal Institution Christmas Lectures which Professor 
Hartridge gave to boys and girls in 1946, and the spirit of entertainment and instruction 
always associated with these lectures is well reflected in its pages. It is essentially an 
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expcrimental book in which a great many phenomena are mentioned. These include 
the spectrum, spectral absorption, colour mixture, interference, diffraction, applications of 


colour in science and medicine, dichroism, the effect of intermittent illumination and so on, | 


Anyone invited to give a lecture on colour will find a wealth of demonstrations from which to 
draw, many of the experiments being illustrated in colour. ‘The theoretical explanations of 
the phenomena do not go very deeply into the subject, but that is hardly to be expected in a 
book of this type. Quite reasonably, Professor Hartridge has not considered this the 
appropriate opportunity to enlarge on his theoretical ideas about colour vision. 

W. D. WRIGHT. 


Cathode-Ray Tube Traces, by HtLary Moss. Pp. 66. 1st Edition. (London: 
Electronic Engineering, 1949.) 10s. 


This monograph, based on articles published by the author in the journal Electronic 
Engineering, should be of interest to all users of the cathode-ray tube. 

The subject matter in general is well presented and the mathematical arguments are not 
undu! comp.icated. 

The first chapter merits separate mention, in that it deals in detail with the various 
forms of Lissajous Figure and, in addition to the basic theoretical considerations, points out 
various practical pitfalls. 

Chapter II is disappointing, both in title (‘ Straight Line Time-Bases ’) and in content. 
It is surprising that the linear sawtooth waveform should be dismissed in a mere section 


of this chapter when various other time-bases of comparatively rare occurrence, such as the 


circular and spiral forms, should merit the whole of Chapter ITI. 

On reading these chapters one gathers, perhaps wrongly, that Dr. Moss has over- 
concentrated his attention on time-bases of sinusoidal origin rather than on those of 
discontinuous form. 

The remainder of the book is mainly devoted to the interpretation of various waveforms 
more usually encountered as the phenomenon under investigation. 

A résumé is given of Fourier series and application is made to several of the common 
pulse waveforms. Beats and the amplitude modulated envelope are also treated. 

The author deserves congratulations for his illustrations, but the bibliography is sketchy 
and would only serve as a first guide to further reading. A. V. LORD. 


Photoelasticity: Principles and Methods, by H. 'T. Jessop and F. C. Harris. 


Pp. vilit+184. 1st Edition. (Lordon: Cleaver-Hume Press Ltd., 1949). 
28s. 


A book on photoelasticity emanating from the Laboratory whence came the classic 
treatise on the subject by Coker and Filon commands immediate interest and attention. 
Work on photoelasticity continues at University College, London, and the senior author 
was the first Chairman of the recently formed Stress Analysis Group of the Institute of 
Physics. 

Authors are not responsible for their publisher’s advertisements, but the assertion on the 
dust jacket that it is a reference work is supported by the authors’ statement in the preface 
that they have given in the text the sources of any original work quoted. Unfortunately, 
and this is the most serious criticism of the book, it is not true. | Statements are often made 
without the name of the worker responsible being given, seldom if ever is a reference 


accompanied by chapter and verse, and there is not even a bibliography from which one 


might hope to guess at the origin. The tabulated data in the appendix would inspire more 
confidence if reference were made to the sources. 

A laboratory such as that of the authors which is the centre of so much work on the 
subject provides an embarras de richesses and the authors seem to have had difficulty in 
selecting and arranging their material. The book gives the impression of having been 
hurriedly put together, and not enough care has been spent in critically examining it before 
printing. For example on page 39 the authors state that a convenient convention of signs 
is the Cartesian; yet on page 43, they give the distances from a principal point of a lens 
system to the corresponding principal focus, H,F, or HF, as both equal to fy; whereas in 
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a Cartesian system one must be the negative of the other. Again on page 63, the definitions 
of the stress-optical coefficients leave the reader guessing as to which refers to the ordinary 
and which to the extraordinary ray, and what is the sign of the difference between them. 
Again on page 60 a compensator is said to be between ‘ crossed ’ polarizer and analyser, yet 
it is not until page 65 that a definition of this term is given. When much is left vague, the 
exact analysis of the mode of working of quarter-wave plates given in the appendix, seems a 
little out of place, especially in a book which sets out to be a guide to ‘‘ those who wish to 
undertake practical work on stress-exploration ’’. 

The book is well printed and the diagrams are clear, but the half tone illustrations are 
poorly reproduced; it is a pity these were not grouped together and printed separately on 
glazed paper; justice would then have been done to the skill of Mr. Storrer the photographer. 
It would seem that the worked examples are invaluable, and the descriptions of the practical 
requirements of the optical bench are admirable, although emphasis should have been placed 
on the need for examination in parallel light, more particularly when the specimens are 
thick and highly stressed. 

It is to be hoped that when the authors prepare a second edition they will expand 
Chapters IV and V on polarization and photoelasticity, and will give references throughout, 
so as to enable those who may wish to go more deeply into the subject to do so with less 
difficulty. A. M. TAYLOR. 


Symposium on Electronics, edited by A. G. Peacock. Pp. xiii+199. 1st Edition. 
(London: Chapman and Hall in collaboration with S.I.M.A., 1949.) 16s. 


The first of a series of annual Symposia on Electronics in Research and Industry, 
organized by the Scientific Instruments Manufacturers Association, was held in London in 
1948. ‘The papers presented at this Symposium are published here in book form, and cover 
a valiety of applications of electronic methods to the solution of instrument problems. 
Recent advances in the measurement of frequency, small displacements, ionizing radiations, 
low gas pressures and sound are described, and accounts are also given of the use of electronic 
techniques in computing, the radiosonde, ultrasonics, spectroscopy, picture telegraphy and 
metal detection. Because of the large number of subjects covered the treatment of each is 
brief and descriptive. For this reason the book is primarily of interest only to those 
non-specialists who wish to have a general summary of a few of the many ways in which 
electronic equipment is being applied, both in research and in industrial production. 
Accounts are given of the discussions which followed the presentation of each paper, but 
the space so occupied could, in the reviewer’s opinion, have been more usefully allocated to 
the expansion of some of the papers. J. M. MEEK. 


Ultra-Violet and Daylight Rays, by J. R. ASHworTH. Pp. x+73. Ist Edition. 
(Liverpool: University Press, 1949.) 6s. 


This monograph presents the results of Dr. Ashworth’s measurements of ultra-violet 
and visible light received at ground level over the 16 years from 1932 to 1947. He has used 
a simple type of photographic measuring technique, which is described and compared with 
other techniques. The results are given in the form of tables and graphs, showing the 
yariation of the average intensities of visible and of ultra-violet light, with wavelength near 
3600 A. measured in arbitrary units. Evidence is given that these very simple measure- 
ments have some real significance, and are not unduly affected by the characteristics of the 
particular place where they were made (Rochdale). 

The results are interpreted as showing that the visible light emitted from the sun is 
roughly constant throughout the sunspot cycle, but that the ultra-violet is more intense at 
times of sunspot maximum. ‘These conclusions are, however, not directly deducible from 
the observations because, according to the author’s theory, the variable ionization produced 
by the ultra-violet in the upper atmosphere in turn produces variable absorption in both the 
ultra-violet and the visible light reaching the earth. 
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It is unfortunate that the theory is left in a very speculative and indefinite form. It is 
suggested that the absorption occurs in Region E of the ionosphere, and is due somehow to 
the ionization of this region; it would be interesting to know in more detail what Dr. Ash- | 
worth believes this mechanism of absorption to be. Reference to his published papers fails. 
to provide any more detail. J. A. RATCLIFFE. 


Les ondes électromagnétiques centimétriques, par LOUIS DE BRoGLIE. Pp. 274. 
(Paris: Editions de la ‘Revue d’Optique Théorique et Instrumentale ’, 


1948.) 800 fr. 


Each year since 1944 M. Louis de Broglie has organized a conference to summarize the 
state of knowledge in some important branch of physics. In 1947 one of these conferences, 
occupying five days, was held to discuss what was known of centimetric electromagnetic 
waves, and the proceedings are recorded in the volume here reviewed. 

The object was not to repeat the numerous accounts which have already appeared of the 
workings of radar, and of the magnetron, which led to the possibility of centimetre radar, 
but to choose nine separate and less well-known aspects of the subject and to ask nine 
separate authors to report on these. These reports are collected in this volume. They 
deal with the subjects of : the exchange of energy in electron tubes; the effects of space 
charge in velocity-modulated tubes; the effects of curvature in waveguides; excitation of 
and coupling to cavity resonators; measurements at centimetre wavelengths; aerials; 
propagation; radiation from the sun and the galaxy; the use of centimetre waves in nuclear 
physics. 

In a conference of this kind it was inevitable that the treatment should differ from subject 
to subject. Some of the papers describe some special point rather fully; those on space 
charge and on curved wave guides are of this type. Others, for example those on aerials 
and propagation, attempt to provide a broad survey without much detail. 

All the subjects dealt with are advancing so rapidly that a report of a conference held 
three years ago must inevitably seem somewhat out of date. Since the study of extra- 
terrestrial radiation was only started seriously after the end of the war, it is perhaps only 
natural that the corresponding section of the report now seems most out of date. It is also 
a little surprising to see this subject, whose advance has depended so largely on measurements 
of metre waves, so fully discussed in a volume dealing with centimetre waves. 

J. A. RATCLIFFE. 


Recent Advances in the Physiology of Vision, by H. Harrripce. Pp. xii+401. 
(London: J. and A. Churchill Ltd., 1950.) 25s. 


Much of Professor Hartridge’s book is concerned with topics belonging to physiological 
optics and colour physics, in which physiologists and physicists have been equally active. 
Eye motions and the dioptric properties of the eye, visual resolving power, spectral lumino- 
sity functions, the directional sensitivity of the retina, quantum notions in vision, colour 
discrimination, the 1: ichromatic and other colour theories, etc., are treated, along with the 
topography of the nerve fibres in optic nerve, optic tract and visual cortex, the electrical 
response of the retina, the effect of light stimuli on the Berger rhythms, the visual pigments, 
electric and magnetic phosphenes, and so on. The author’s object has been to survey a very 
wide field, giving enough of the older work to make clear where advances have been made. 

It is probable that no worker in the subject of vision will read this book without wanting 
to take issue with Professor Hartridge on some of his contentions. But a work of this kind, 
particularly one dealing with vision, can scarcely be criticized for not presenting a coherent 
and generally accepted picture. It can be judged successful if it draws the reader’s attention 
to new results and gives some idea of how people are interpreting or speculating about them. 
Professor Hartridge’s book does all these things, and is, perhaps, particularly strong on the 
‘ speculating’. The style is easy, but the arrangement of material leaves something to be 
desired. In a number of cases the same subject is dealt with in widely separated sections. 
For example, eye-motions and their measurement are discussed in Section 8 (p. 64) and in 
Section 54 (p. 329), but there seems no justification for dividing the material in this way. 
Among several topics omitted which might have been expected to earn a place may be 
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mentioned particularly W. J. Schmidt’s polarized-light analysis of the structure of the outer 
segments of the retinal rods, an analysis largely confirmed by the electron microscope 
photographs of rods recently obtained by Sjéstrand. 

Of the author’s handling of visual theories it is sufficient to say that most physicists would 
demand more precise and quantitative reasoning to support some of the statements made. 
A few comments on other aspects of the presentation may be added. Although the crux of 
Hecht, Shlaer and Pirenne’s work on quanta and visual thresholds is the relation between the 
“sharpness ’ of the threshold and the minimum number of quanta involved, no mention is 
made of this relation in the two pages devoted to the work in question. The results of 
Mandelbaum and Mintz on the dark adaptation of the colour receptors, briefly described on 
p. 50, are stated to refer to rod vision. This is clearly a mistake, which is aggravated by a 
suggested connection of the results with dark-adaptation measurements by Miles which 
certainly do refer to rod vision. In deploring the lack ot information on the variation of 
visual acuity with pupil size the author seems to have overlooked the extensive measurements 
of Fabry and Arnulf published in 1937 (La vision dans les instruments, Editions de la Revue 
d@Optique). An explanation is given of how the binocular threshold may be lower than the 
monocular threshold purely on account of ‘ probability summation ’ (pp. 237-238). This is 
immediately followed by a statement implying that any such difference must arise from 
physiological summation. Itis asserted (p. 203) that the coefficient curves for protanomalous 
subjects show no negative coefficients in any part of the spectrum when primaries 650, 530, 
and 460 my are used, and the curves shown in Figure 160 are cited in support. The figure 
is obviously wrong, as at wavelengths below 460 my the coefficients do not add up to unity, 
and the statement itself is incorrect. 

A number of other criticisms of the same kind could be made, but taken as a whole the 
book provides a collection of facts and current speculations on vision which will prove of 
value to all students of the subject. iis Ge S, 


Meéthodes de calcul dans des problémes de mécanique. Pp.102. (Paris: Centre 
National de la Recherche Scientifique, 1949.) No price. 


In the spring of 1948, two international conferences on methods of calculation in 
problems of mechanics were organized at Marseilles and Paris by the C.N.R.S.; the 
papers contributed to these conferences have been collated and published in this volume. 
‘The contributors and topics are as follows: J. Valensi (Marseilles), ‘“ The role of applied 
mathematics in engineering’’; D. N. de G. Allen (London), “‘ Relaxation methods and 
problems of frameworks ”’, ‘“‘ Relaxation methods and the solution of differential equations ”’, 
together with ‘‘ Supplementary notes on the application of relaxation methods ’’; Th. Vogel 
(Marseilles), ‘“‘ The escalator method for the calculation of eigenfrequencies ’’; M. Picone 
(Rome), ‘‘ New points of view in harmonic analysis’; L. Couffignal (Paris), “‘ Calculating 
machines for harmonic analysis according to the method of H. and Y. Labrouste’’; J. M. 
Burgers (Delft), ‘‘ Problems related to the theory of turbulence’’; L. Malavard (Paris), 
*“ Some recent applications of the method of electrical analogies’’; A. van Vijngaarden 
(Amsterdam), ‘‘ Potential flow about a solid of revolution”; F. H. van den Dungen 
(Brussels), ‘‘ The application of the calculus of variations in fluid mechanics ”’; L. Couffignal 
“<The role of numerical calculation in scientific and technical research’’. 

The papers naturally vary in style and standard but they are, without exception, lucidly 
written and the book can be thoroughly recommended to a physicist or engineer who wishes 
to acquire some idea of modern techniques used in the realm of classical mechanics. R. M. D. 


Glass-to-Metal Seals, by J. H. Partripce. Pp. xii+238. Ist Edition. 
(Sheffield: Society of Glass Technology, 1949.) 25s. 


This is a comprehensive and authoritative guide to the art of sealing metals to glasses. 
The production of graded glass-to-glass seals is therefore also dealt with in detail. The 
problems of stress distribution, so vital in determining success or failure, are ably and 
adequately discussed. 

This book should be available to all workers, both scientific and technical, in every 
physics and chemistry laboratory; it is first-class. Cork 
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ABSTRACTS FOR SECTION A 


A Magnetic Cosmic-Ray Spectrograph with Counter Recording, by B. D. Hyams, 
M. G. My rot, B. G. OweEN and J. G. Witson. 


ABSTRACT, An instrument is described in which the momentum of single cosmic-ray 
particles is measured, the particles leaving the instrument in a sufficiently collimated beam 
for use in subsequent experiments. The ‘maximum detectable momentum’ of the 
apparatus is about 3 x 10!° ev/c., almost all particles are collected within 10° of the zenith 
and in directions individually defined to +1°, the output is about 400 particles/day, and the 
proportion of spurious records is of the order of, or less than, 0:1%. 


A Magnetic Cosmic-Ray Spectrograph with Counter Recording—I1: The Electronic 
and Recording System, by B. G. OweEN. 


ABSTRACT. The paper describes the electronic selection and recording systems of the 
spectrograph and the performance and overall serviceability of the equipment as shown in 
continuous operation over about one year. 


Observations on Cosmic-Ray Penetrating Showers at High Altitude, Sea Level 
and Below Ground, by E. P. GrorcE and A. C. Jason. 


ABSTRACT. Penetrating showers have been investigated, using Geiger counters, at a 
depth of 30 m. below ground, at sea level and at an altitude of 3,457m. The influence 
of the geometry of the absorbers was studied. The transition curves for local showers, 
corrected for this geometric effect, give the following collision lengths of the generating 
particles; Pb, 180-+40; Al, 85+15; Paraffin ~80 gm/cm?. These values correspond to a 
collision cross section close to the geometric cross section of the nuclei. From the results 
obtained for extensive penetrating showers, it is concluded that most of the penetrating 
particles in these showers appear to be generated in the atmosphere. In 40 days, no pene- 
trating showers were recorded below ground. Using photographic plates, the density 
distribution of shower tracks associated with nuclear disintegrations was found to be similar 
at sea level and 3,457 m. 


Slow Mesons in the Backward Flux of the Cosmic Radiation, by D. M. Ritson. 


ABSTRACT. Using the method of ‘ delayed coincidences’ measurements have been 
made at sea level on the flux of slow mesons travelling backwards from the ground. 


Fluctuations of the Electron—Photon Cascade—Moments of the Distribution, by 
L. JANossy and H. MesseEL. 


ABSTRACT. 'The second moments of the electron—photon cascades are evaluated. It is 
found that the mean square fluctuation at the cascade maximum is that expected for a 
Poisson distribution or less; before and after the maximum the fluctuation exceeds greatly 
that expected for a Poisson distribution. The correlation coefficient between electrons and 
photons is zero for small thicknesses or near the cascade maximum. For other regions 
the correlation coefficient is positive. Both mean square fluctuation and correlation between 
electrons and photons can be understood qualitatively in terms of fluctuation of effective 


depth. 


The Effect of Nuclear Structure on the Elastic Scattering of Fast Electrons, 
by L. R. B. ELTON. 


ABSTRACT. The cross section for the scattering of electrons by atomic nuclei 1s 
investigated at energies for which the nuclei can no longer be treated as point charges. 
Two simple nuclear models are used. A general expression is obtained using Born’s 
approximation, and an exact numerical calculation is carried through for 20 Mey. electrons 
scattered by gold nuclei. It is concluded that at this energy considerable deviations from 
the formulae which treat the nucleus as a point charge should be expected, and that these 
should furnish some information about the charge distribution within the nucleus. 
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Non-physical Solutions in Classical Finite Electron Theory, by J. IRVING. 


ABSTRACT. An approximate linear form of the Peierls—McManus equations of motion | 
for an electron is shown to yield runaway solutions for a particular influence function, — 
which otherwise satisfies all the conditions imposed by the theory. The general problem — 
is to discuss the existence of complex roots of a transcendental equation for an arbitrary 
form of the influence function. It is shown that this equation always has roots. ‘These, | 
however, may be such that the corresponding motion of the free electron is of a damped ~ 
nature. It has not been found possible to construct a function which avoids runaway | 
solutions or to prove whether or not such a function exists when all conditions are satisfied, | 


The Band-Spectrum of Barium Oxide, by A. Lacerovist, E. Linp and 
R. F. Barrow. 


ABSTRACT. The band-spectrum of BaO between 5000 and 7000 a., which consists of 
a single 1X—1X system, has been the subject of a new rotational analysis in which the follow- 
ing bands have been examined: 5,0, 4,0, 3,0, 3,1, 2,0, 2,1, 1,1, 1,2, 0,2, 0,3 and 0,4. 
‘The constants derived are : 


B,’ =0-3124,—0-0013,(v” +3)—0-0000,(v"+3)2 BB,’ =0-2584—0-0011,(v’ + }) 


De=26-5 10m DI=28>105* 
re =1-940 x 10-8 cm. re =2:133 x 10~* cm. 
We =669°8,, Me Gs SOS: we ~500, x”w ‘=1", 


Vo,o= 16,722-2; cm}. 
In the upper state some 16 perturbations have been found. They appear to arise from 
interactions with at least four electronic states or substates. Fairly complete information 
about one of the perturbing levels has been obtained. ‘The constants of this state are: 
B,=0-2254—0-0013,(v+4), 
Gy=450-4(0+4)— 2+ D2, 
Vo,90=17476-, em. 


The band-constants of the other perturbing states are very similar to those given above. 


Band-Spectrum of ZrO, by M. Arar. 


ABSTRACT. The spectrum of ZrO has been photographed in the ultra-violet, the 
visible and the infra-red, from arcs running under special conditions. This has made 
possible the identification of three new systems in the ultra-violet, as well as two less 
conspicuous ones in the infra-red. It has also made it possible to extend the blue system « of 


ZrO. 
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